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COMOC II: TWO-DIMENSIONAL AERODYNAMICS SEQUENCE, 

COMPUTER PROGRAM USER'S GUIDE 
By 

P.D. Manhardt, J.A. Orzechowski 
And 

A.J. Baker 

COMPUTATIONAL MECHANICS CONSULTANTS 
Knoxville, TN 


SUMMARY 


The COMOC finite element fluid mechanics computer program system is 
applicable to diverse problem classes. The two-dimensional aerodynamics 
sequence has been established for solution of the potential and/or viscous 
and turbulent flowfields associated with subsonic flight of elementary 
two-dimensional isolated airfoils. The sequence consists of three 
specific flowfield options in COMOC for two-dimensional flows. These 
include the potential flow, option (2DPF), the boundary layer option (2DBL), 
and the parabolic Navier-Stokes option (2DPNS). By sequencing through 
these options, it is possible to computationally construct a weak-inter- 
action model of the aerodynamic flowfield. This report is the user's 
guide to operation of COMOC for the aerodynamics sequence. 


INTRODUCTION 

The COMOC finite element computer program can solve various differ- 
ential equation descriptions for fluid flows of practical interest. The 
two-dimensional aerodynamics sequence has been established to allow deter- 
mination of the inviscid and irrotational flow about an isolated airfoil, 
and the resultant pressure distribution including the trailing edge wake 
trajectory. The potential flow shape may be the actual airfoil, or an 
"effective" inviscid shape established by inclusion of viscous and turbulent 
flow effects in the potential solution. These corrections are determined 
by solution of the boundary layer equations for the viscous attached flow, 
and the parabolic Navier-Stokes equations for the merging of upper and lower 
surface turbulent boundary layers into the trailing edge wake. Under the 
assumption of weak-interaction, the viscous flow computations are executed 
within the axial pressure (gradient) field established from the inviscid 



flow solution. The program outputs all viscous flow parameters of interest. 
The computed distribution of displacement thickness can be employed to 
augment the inviscid airfoil shape, and a new potential flow pressure dis- 
tribution then established. Therefore, by cycling through the various flow- 
field options in COMOC, it is possible to build a weak-interaction solution 
for the aerodynamic shape of interest. 

The two-dimensional aerodynamics sequence in COMOC is constructed upon 
three distinct program options. Each option is basically a computational 
sequence that allows solution of the differential equation (s) governing 
specific classes of fluid flow. The potential flow option (2DPF) establishes 
the distribution of velocity perturbation function for inviscid irrotational 
flow about an isolated two-dimensional airfoil. It can be instructed to 
call the pressure coefficient routine which computes Cp from the potential 
distribution on the airfoil and on the trailing edge wake. The boundary 
layer flow option (2DBL) established solution for turbulent two-dimensional 
attached flow in a given Cp distribution. It can run two (opposed) boundary 
layer solutions with distinct Cp as well, as occurs with upper and lower 
surface flows on the airfoil. In both cases, closure for turbulence can be 
accomplished using mixing length theory (MLT), or a second order model built 
on turbulent kinetic energy (TKE) and solution of additional differential 
equations. Upon special input specification, an economical integral boundary 
layer solution can be computed as well, for the entire airfoil including 
transition from laminar to turbulent. 

The merging of the upper and lower surface boundary layers at the trail- 
ing edge into a turbulent downstream wake is an important part of the inter- 
action flowfield. Solution for the wake flow is accomplished using the para- 
bolic Navier-Stokes (2DPNS) option in COMOC. The lateral boundaries for 
this flow are the extended upper and lower surface inviscid flows and the 
solution is marched down the wake. Closure for turbulence must now be the 
TKE model; initializing levels are obtained from the boundary layer solutions. 
Both the 2DBL and 2DPNS options provide detailed downstream distributions of 
important integral parameters, e.g., skin friction (if appropriate), dis- 
placement thickness, shape factor, etc. In the wake region, these parameters 
are referenced to the downstream projection of the trailing edge included 
angle bisector. The displacement thickness is usually interpreted to define 
the effective boundary of an inviscid flow. Therefore, the computed distri- 
bution of displacement thickness on the airfoil surface and wake can be em- 
ployed to define a new effective inviscid shape, and the sequential solution 
procedure repeated. This constitutes the weak-interaction solution algorithm, 
which the aerodynamics sequence in COMOC is capable of executing. 

COMOC is undergoing a rapid rate of growth and the features of the pro- 
gram discussed herein are newly developed. The theoretical foundation for 
the finite element solution algorithm, applied to the components of the weak 
interaction solution, and detailed analysis of accuracy, convergence, and 
verified capabilities are presented in reference 1. This report is intended 
to guide the program user in correct problem specification and data deck 
preparation. It contains a discussion of the organization of COMOC, including 
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a complete description of all real and integer arrays and all FORTRAN varia- 
bles. Brief mathematical descriptions of^the problem statements are given 
in the next section, including boundary condition constraints. An overview 
of the finite element algorithm form is presented. The next section dis- 
cusses data deck preparation, including detailed output for comparison 
standard test cases. The final section presents the organizational and 
logical structure for the program. 


PROBLEM DEFINITION 


Aerodynamic Potential Flow 

The 2DPF option in COMOC solves the elliptic boundary value statement 
for the perturbation potential function, which is the Laplacian, i.e.. 


L(4>) = 




3x^. 3x^ 


0 


( 1 ) 


The definition for (|) , in terms of the reference freestream velocity at 

angle of attack and the local velocity vector Uj is 



Equation (2) provides the flow tangency boundary condition used for solution 
of equation (1). Denoting fi-; as the unit normal vector defining any 
impervious surface, the boundary condition is 


d(f) ^ ^ ✓v 


(3) 


In the freestream, sufficiently remote from the airfoil surface, the flow 
is undisturbed. Hence, equation (3) vanishing identically can account for 
this as well as the Kutta condition to admit angle of attack. 


Turbulent Boundary Layer Flow 

COMOC can execute this option using a coordinate stretching transfor- 
mation to economically account for boundary layer growth (see ref. 1, 
equation 94). Denoting as the time-averaged mean flow velocity vector, 
and assuming constant density, the 2DBL option solves the differential 
equation system 
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(7) 


In equations (4)-(7), v® is the effective kinematic viscosity defined as 

V® = Re~ + V ^ 

(8) 

Here, Re is flow Reynolds number and is the turbulent eddy viscosity, 
defined as 


’ 

au. 


1 

< 

8X2 

(Cvk^e-i 


(MLT) 

(TKE) 


(9) 

( 10 ) 


In the MLT specification, u is the VanDriest damping and i is the mixing 
length. 

Boundary conditions are requiredfor all four variables, equations 
(4)-(7), at the wall and all except U 2 at freestream. The freestream 
condition is simply vanishing gradient. The no-sTip constraint is used for 
ui and U 2 may have a specified value at the wall. The TKE closure is 
invalid at a no-slip wall, since the local flow is not fully turbulent. 

COMOC employs a sub-layer model to internally evaluate appropriate boundary 
values for both k and e , at a user-specified value of y+ (see ref. 1). 


N 
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Turbulent Wake Flow 


The 2DPNS option in COMOC can be entered directly under a restart from 
the 2DBL option. The 2DPNS equation system is. 



L(pu^.) 

L(k) 


L(e) 


3u, 


pu 


j 3X . 


+ ^ - 


3Xi 


3x, 


r8u,. ^ suj' 


3x, 




= 0 




V 


3x^ 3x^ 


^UjC) - 


3 

V® 3k 


P>^k 

e 

= 0 

3 

V® 3e 

axj 

Pr^ 3Xj^ 

lUj 

(Xj, 

+ c|e^ 


= 0 


( 11 ) 

( 12 ) 


(13) 


K14) 


where 1 < i , j <2 and i = 2 only. The second of equations (12) is now 
solved for > with equation (11) providing the freestream boundary 
conditions in'^terms of Cp, top and bottom. Since the wake is an unbounded 
flow, the slip surface has vanished and the solution domain for all variables 
is closed by vanishing freestream normal gradient. Further, the TKE equation 
system is everywhere valid. 


Finite Element Solution Alqbrithm 


The appropriate differential equation systems are established. Each is 
a special case of the general, second order elliptic boundary value specifi- 
cation 



3x, 


K-^ 


+ f. 




+ fr 


In , ,-^1 

l’3x. 


= 0 


(15) 
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where q is interpreted as a generalized dependent variable, i.e., 
q = Uj,.U 2 ,k,e }. The boundary condition statement applicable for all 
numbers or the set q is 

£(q) = a^^^q + + a^^^ " ° 1 (16 

The finite element algorithm for equations (15)-(16) is established 
using the method of weighted residuals. It is 



{(t)(x£)} L(q*)dx - X 


{<|)(Xj^)}il(q*)dT| 





(17) 


where S is the finite element assembly operator, and within the finite 
element'^domain Rf^ , the finite element approximation qjj^ to the solution 
q is 


(x^) H {f(x.^)}T {Q(x)};^ 


(18) 


where {.f(Xjj;)} are polynomials written on coordinates spanning Rjp* Upon 

application of a Green-Gauss theorem, the globally assembled algorithm for 
the representative differential equation system (15) -(16) is 



In equation (19), the a^^^ are coefficients in the boundary condition 
statements. Of particular interest, note that for q = (f) , 


a^^^i h. £. (20) 

which introduces the non-homogeneous boundary condition constraint directly 
into the algorithm. 
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DATA DECK PREPARATION 


The input facilities in the COMOC weak interaction viscous-inviscid 
interaction flow program are highly sophisticated and greatly simplify data 
deck preparation and modification. The program sequentially scans the data 
deck and operates on command data cards as they are encountered. Numerical 
data required for each cotmand operation is input in free format on cards 
directly following the command card. Command operations can cause vectors 
to be filled, initiate a series of solution operations or specify output 
formats and titles. Command card sequence is quite flexible and care has 
been taken to ensure that most operations which must be performed sequen- 
tially are specifiable under one command name. Guidelines for numerical 
data preparation^command card sequencing and utilization of the output are 
described in this section. Copies of actual data decks used for solution 
of the test cases are given in the appendices for reference. 


Structure and Guidelines 


The COMOC data deck for weak interaction flow solution is divided into 
six sections for description. Exclusive of machine related job control 
cards, the six sections consist of a Fortran main routine and accompanying 
subroutines O'f any), namelist data, geometric description, output format 
specification, boundary and initial condition data and solution command 
cards (Fig. 1). Each of these data and its subset is preceded by a 
command card which directs a program activity which when completed returns 
control to the next data card. The program operates in a dynamic storage 
mode and the function of Main is to allocate sufficient storage for the IZ 
array which is internally sized as a function of the number of finite elements 
requested for a specific problem. The namelist section of the deck is used to 
specify scalar integer and floating point data utilizing the Fortran namelist 
option. The data is read in namelists NAMEOl and NAME02, respectively, and 
stored in the arrays lARRAY and RARRAY. Execution of the namelist read is 
initiated via commands FENAME or NAMELIST. 

The geometric description section contains the data required to generate 
a finite element grid suitable for solution. Command FEDIMN dynamically 
dimensions the arrays required by the analysis according to the number of 
nodes (NODE) specified in namelist NAMEOl. The finite elements are generated 
via command card LINK2 (14) for viscous flow and LINKl (09) for potential flow. 
Finite element numerical data follows these cards as noted in the data deck 
description. 


COMOC employs a highly adaptive output routine which allows for data 
specification of the scalar and array variables to be printed, scale factors 
to be applied to each variable and titling information to head each variable 
list. Each is specified under a command name in the output section (Section 
IV, Fig. 1) of the data deck. The program operates in non-dimensional units 
and the data specified scale factors are utilized to present the results in 
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Solution Procedure 




a consistent set of units. A reference length parameter (REFL in NAME02) 
is also available to scale all output to a problem reference length, i.e., 
airfoil chord length. In addition, command names are available for speci- 
fication of problem identifying titles to be printed at various strategic 
locations such as the beginning of each set of printed output. 

The fifth section of the data deck contains the required boundary and 
initial condition data which are specified at the solution nodes. The finite 
element method easily handles mixed boundary conditions, hence, both fixed 
and Neumann type are allowed. Parameter tables such as C vs x may also 
be specified in this section. Initial conditions are ^ required only 
for the viscous flow calculations. The inviscid (Cp) solution boundary 
conditions are automatically generated in subroutine ARFOIL , hence, need 
not be specified. 

Section VI contains the commands for the solution proper. For poten- 
tial flow, the Laplacian elliptic boundary value problem is solved on the 
perturbation potential function. Pressure coefficient Cp is subsequently 
evaluated over the airfoil and in the wake region by differentiating ({). 

For viscous flow solutions, the generated Cp curve is differentiated in 
DPDXTB . Upon command QKNINT in the solution section, integration of the 
2DBL equation system over the airfoil or integration of 2DPNS in the wake 
region is performed subject to the applied pressure gradient. 

A complete listing of the allowed command names and a brief description 
of their use is given in Table 1. During the progression from sections one 
to five in the data deck, various execution and data management commands are 
required in addition to numerical data specifications. These are more fully 
explained for particular sample test cases which follow. 


Free Format 

Most of the numerical data other than command data specified in the 
COMOC input deck may be input in free format. Data delimiters may be blanks 
or commas. The end of a data set is indicated by a T or blank card for 
numerical data and a DONE beginning in card column 1 for literal data. 
Exceptions to the rule are namelist data which utilizes the standard Fortran 
Namelist option and certain special card types which combine literal and 
numerical data. Command cards are an example of this type and the restric- 
tions imposed are that the command name begin in column 1 and numerical data 
begins beyond column 10. 

Several features which greatly simplify sequential and repetitive data 
specifications are available in free format. For example: 
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Repetitive Numbers 


12. 5*7. 


Fills Array 

Repetitive Sequence 
(one per card only) 
Fills Array 

Skip P locations 
Fills Array 

Increment by a constant 
Fills Array 
Exponential Notation 
Fills Array. 


12. 7. 7. 7. 7. 7. 
2(5. 2. 4. 

5. 2. 4. 5. 2. 4 
10. 12. 3*P 22. T 

10. 12. A A A 22. 

5* 50 10 T 
10 60 110 160 210 
6. lO.E-2 14.E-4 T 
6. .1 .0014 


Restart 

The airfoil viscous solution proceeds in marching fashion along the 
upper and lower airfoil surfaces toward the trailing edge and into the wake 
region. If the solution fails to reach the final station in the allotted 
Cpu time, it is desirable to restart the solution where it ended rather than 
begin again. The RESTART and SAVETAPE commands permit this by writing the 
entire stored array on a user specified file at each requested print interval. 


Sample Test Cases 

Three test cases were selected for checkout of the COMOC solution options 
which are coupled under a weak interaction assumption. Test case results are 
presented in the accompanying theoretical report (ref. 1). Potential flow and 
pressure coefficient (Cp) were calculated for flow over a NACA-0015 airfoil 
with modified leading edge region as noted in Ref. 2. Secondly, a viscous 
boundary layer solution was obtained for the geometry of Bradshaw as reported 
in the Stanford Conference Proceedings (ref. 3). The wake region was analyzed 
for a Joukowski 12% thick airfoil at 6° angle of attack. Data deck setup for 
these cases is described in detail in the following sections and a copy of the 
complete data deck for each is presented in the appendices. 
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TABLE 1 


Command Name Input ^ 

f3:.' 


A CCNTFQL CARO W[TH TH6 FGLLCWING FARAMETERS IS READ IN - 


PAR AMETE P 

FORMAT 

CARD 

CCLS 

CESCR IPTION 

VI 

A8 

1 - 

8 

CONTROL 

VARIABLE. 

AMUL 

FREE 

AFTER 

8 

NMULI 1 ) 

= NX 





NMULI 2) 

= NPCD 





AMUL 13) 

= NREPET 





AMULI 4) 

= NR TAPE 




TG 

AMUL(N) 



IF KOUNF = I IN NAMEOlt TFEN TFE ENTIRE INPUT CARO IS PRINTED 
IMMECIATELY AFTER BEING READ FCLLCNEC RY THE DATA THAT IS BEING 
STORED Al.CNG WITH TFE DATA'S ENTRY POSITION IN THE \l ARRAY 
(SEE IPLACE AND AFLACE) 

THIS ROUTINE LCCKS FIFST FOR A MATCH CF VI WITH CERTAIN KEY WORDS 
WHICH WII,L EITHER CALSE A SUBROUTINE TO EE CALLED OR PROGRAM 
FLCW TO CCCUR. 


THE KEY WCPCS THAT ARE SC ANN EC ARE - 


ABSVAL Nl A2 N3- 
l BLANK) 

CCVCC 
COMTITLE 
DESCPIFT NX 


END 

EXIT 

2DPF 

3CRR 

3CPNS 

2DNS 

FEDIMN 

F ENA ME 


SET RZ(IZ(M2)) ’ ARSIRZI I2IN3) ), I = I, Nl 
RETURN TO SCAM ANOTHER CARD. 

PRINT THE COMGC TITLE PAGE TWO TIMES. 

READ A TITLE CARD WHICH WILL APPEAR ON CGMOC . 

CALL CSCRPT AND PROCESS ACCORDING TO NX. 

NX 

FLANK - READ AND WRITE INFQRjMATION CARDS. 

203 - READ TITLES FOR OEP. VAR. OUTPUT HEADINGS. 
20A - READ DESCRIPTIVE TITLE FOR HEADING 
AT BEGlNi'lINO OF OUTPUT HEADER. 

332 - READ TITLES FOR PARAMETERS PRINTED IN 
THE OLTPLT HEADER. 

RETURN CONTROL TO MAIN PROGRAM, RESET ARRAYS 
AND RETURN TC I301NFT. 

CALL EXIT. 

set IOL = 0 AND IPHI = 1. 

SET FBI. = 1. 

SET IBL = I AND N3DPNS = 1. 

SET IBL = 0, 

call DIMEMSIONALIZATION ROUTINE FEDIMN. •- 
CALL FENAME TO SET DEFAULT SCALARS AND THEN 
CALL NMELST TO READ IN NAHEOl AND NAME02 
namelists. 


lARRAY Nl, N2, A3, AA, ETC. 

SET lARPAY(Al) = N2, IARRAYIN3) = N4 , ETC. 
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TABLE 1 Contd 


ICCND 


- 

CALL ICCND TC PRINT REAL AND INTEGER SCALARS. 

INPUT Nl 


- 

SET INPUT UNIT TO Nl. 

INTEGER 



ALLCWS NFW VALUES TO READ INTO A SEQUENCE 
CF LOCATIONS IN THE BORDER, I PLACE AND NLOC 
VECTORS. 

KBNO 

MX 

- 

ENTER FIXED NCOES FCR OEP. VARIABLE NX. 




CALL GETBMO 

KBAO 

NX 

1 - 

ENTER FIXED NODES ANC/CR BOUNDARY 
CCNOITICNS FOR OEP. VARIABLE MX. 




CALL GFTfiCD 

LINK 1 

NX 

- 

CALL LINKl(N'X) 

LIAK2 

NX 

- 

CALL LINK2INX) 

LINK3 

NX 

- 

CALL L INK3 (NX ,DUPRYl ,DUMMY2 ) 

L INK4 

NX 

- 

CALL LINKAINX.K) 

LINKS 

NX 

- 

CALL LINK5(NX) 

MATSLM 

N 1 , 

iM2, N 

3 , N4, NS 



- 

CALL MATSUM ( RZIIZ(N2)), RZIIZ(N3)), RARRAY(N4), 




RZnZ(NS)) , Nl ) 

NAMEL 1ST 



CALL NMELST TO READ IN NAMEOl AMD NAME02 
NAVEL I STS. 

POUMP Nl 

N2 

N3 - 

CALL POUMP ( I2IIZ(N2)), IZ(IZ(N3)), Nl ) 

PLUS M 

N2 K3 . . 

. ETC. 

READ Al 

A2 

N3 A4 

READ ( Ml ) RZ( IZ(M3)+M4) ,1=1, N2 

SAVFTAPE 

Nl 

N2 - 

SAVE CUTPUT CN UNIT Ml, REWIND AFTER 


MOD ( KOUMT, N2 ) =0. 

SETVAL Nl N2 N3 A4 A5 

- CALL SETVAL ( R2(IZ(N'2)), RZ(IZ(M3)), 
RARP.AY(M4), RARRAYINb), Nl ) 

SORT M N2 ... - R/!RRAY(Ml) = S CRT < R ARRAY ( M2 > 

RARRAY(.M31 = SCR T ( R ARR A Y ( M4 ) ) 

... FTC. 

SQRT Nl N2 M3 - R Z ( I Z ( N2 ) = SGPT ( R Z ( I Z ( N3 ) ) > , I = 1, Nl 
VYYENO NX - CcNOTFS £N0 OF INPUT FOR OEP. VAR. NX. 

WRITE M N2 A3 A4 

WRITE ( Ml ) RZ( IZIN3)+N4) , I = 1 , N2 


IF ACNE OF THE ABOVE SITUATIGMS OCCUR, 

THE VECTORS BCRCER AAC VALUE ARE SCANNEC UNTIL A MATCH IS FOUND 
ANC the LOCATION IS STORED IN THE PARAMETER 'K'. 
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TABLE 1 Contd 


OESCPIPTIOM OF INTEGER TTPS INPLT 


EOPCER IS A VECTOR OF CCNTPCL NA.«ES hHICH IS SCANNED WITH THE USER 
INPUT CCNTPOL FOP INTEGER INPUT. 


I Z ( I PLACE IK) ) = LCCATICN IN TFE IZ ARRAY AT WHICH TO BEGIN 
STORING INTEGER ENTRIES. 

I ARRAY INLCC (K) ) = AUKEEP CF ENTRIES STOREC STARTING AT 

TZ( IPLACE(K)). 

IF NX .NE. -I, CALL CET3ND TO ENTER INTEGER DATA. 

IF NX .EC. -1, CALL AOCCEL TO ENTER INTEGER DATA. 

SEE GETBND FOR INFLT DESCRIPTION. 


1 

eCROER 

IPLACE (K 1 

NLOC (K ) 

PPFINITION 

3 

THICK 

68 

93 

ELEMENT THICKNESS VECTOR. 

7 

r P I N T 

5 

3 1 

SCLLTICN SEQUENCE VECTOR. 

12 

L INKC ALL 

121 

125 

LINK NCS. TO BE CALLED AT END OF QKNUI N, 

13 

SPECIES 

31 

121 

VARIABLE NOS. FOR SPECIES TO BE RUN. 

U 

lORULT 

123 

67 

OLTPLT VARIABLE MULTIPLIER FROM RARRAY. 

15 

IOSAV6 

124 

60 

VARIABLE LIST TO BE DISPLAYED AT OUTPUT. 

16 

CNTPTS 

127 

47 

CCNTCUP NODES TO BE USED IN 
CONTES, DFCFBL, TRBTHK, WLFLXS, ETC. 

17 

CNTMDS 

128 

128 

NO. OF NODES IN EACH CONTOUR LINE. 

13 

IBCRC 

38 

131 

COLN^ER-CLOCKWI SE LIST OF BOUNDARY. 

21 

ICNUMB 

131 

142 

LIST OF ENTRIES IN RAP.RAY TO BE 
DISPLAYED AT START OF EACH OUTPUT. 

22 

NPAR A 

135 

67 

LIST OF MULTIPLIERS IN RARRAY USED 
TO MULTIPLY IONUMB ENTRIES. 

26 

NX 

117 

67 

NO. OF SUBDIVISIONS / SUPER ELEMENT 
ALCNG DIRECTION 1. 

27 

NY 

1 18 

67 

NO. OF SUBDIVISIONS / SUPER ELEMENT 
ALCNG DIRECTION 2. 

28 

ELEMENTS 

26 

14 

READ IN ELEMENT NODE CONNECTIONS. 


NOTE - 

NLGC(K) OR NVOC(K) = 67 IMPLIES PRESET LENGTH IS NOT CHANGED. 
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TABLE 1 Contd 


CESCRIPTIGN CF PE«L TYPE IMPUT 


VALUE tS A VFCTOP OF CCMFCL ^'A^•ES F.HICH IS SCANNED WITH THE USER 
INPUT CCMTPOL FOR REAL IMPUT. 


I ZINPLACE (K 1 ) = LCCATICN IN TFE IZ ARRAY AT WHICH TO BEGIN 
STORING PEAL ENTRIES. 

lAPRAYINVOC(K) ) = NUA'eSP GF ENTRIES STORED STARTING AT 

IZINPLACEIK)} . 

ROUTINE REDRSL IS CALLED AT THIS TIRE TO ENTER DATA. 

THESE CCNTRCL CAFCS CAN CCNTAIN A GFCUP OF MULTIPLIERS FOR 
THE ENTERED CATA . 

E. G. 

VYV 3 -lOtT.n -27 

1.2 0.9 3.7 T VORTICITY INPUT 

THE PROGRAM SETS AMLLT = -100.0 * RARRAYI3) / RARRAY127) 
THEN VYY ( U = AMIJLT * 1.2 
VYY( 2) - AEULT * 0.9 
VYYOI = AMULT * 3.7 


K 

VALUE 

NPLACE ( H) 

NVGC(K) 

’ ' — " " — 

fHF INITION 

2 

VU3P0S 

65 

178 

X ST.4TICN FOR VARIABLE GRID CHANGE. 

A 

VU3VAL 

66 

67 

SCALE FACTOR FOR VARIABLE GRID CHANGE. 

5 

VTHI CK 

70 

67 

value CF ELEMENT THICKNESSES. 
DEFAULT = 1.') / ALC 

6 

VRHO 

8A 

67 

DENSITY AT NODE POINTS. 
CE FAULT = R HO INF 

7 

vttab 

19 

5S 

TAfiLF LOOK-UP TEMPERATURES. 
CEFAULT = TOFINF 

a 

vcptab 

IB 

67 

TAELE I.CCK-UP SPECIFIC HEATS. 
CEFAULT = CPOINF 

9 

VXICCR 

B9 

16 

Xl-CCr,RCIiMATES AT NODE POINTS. 

10 

VX2C0R 

90 

16 

X2-C0CRDINATES AT NODE POINTS. 

11 

VH 

79 

67 

ENTHALPY DISTRIBUTION AT NODE POINTS. 
DEFAULT = l.O 

16 

VPRESS 

91 

67 

PRESSURE VALUES AT NODE POINTS. 
DEFAULT = PINE 

17 

VSCHN lOF 

11 A 

67 

SCHMIDT MO. DIST. AT NODE POINTS. 
CEFAULT = SCT 

le 

VYV 

32 

67 

CEFFNOENT VAR. CIST. AT NODE POINTS. 

19 

VTEMF 

85 

67 

TEMPERATURE DIST. AT NODE POINTS. 
FEFAULT = TOFINF 

22 

VTK 

38 

67 

THICKNESS OF ELEMENTS IN THICK VECTOR. 
CEFAULT = 1.0 / ALC 
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TABLE 1 Contd 


23 

VSUTILD 

1 33 

67 

STLDVR, 5TL0TR, STLOCR, STLOEX, STLCON 
EMRIES FOR SUTHERLANDS LAW. 

DEF. .1163E-4, 494.0, 204.0, 1.5, 0.0 

24 

VPRA^CTL 

134 

67 

PRA.VOTL NO. 01 ST. AT NODE POINTS. 
CEFAIJLT = PR 

25 

VX3ST 

139 

161 

DOWASTREAM STATIONS AT WHICH PRESSURE 
IS DEFINED. 

Zt 

VPVSX 

1 40 

67 

PCtASTPEAN PRESSURES AT VX3ST. 

27 

VEPS I LCN 

136 

67 

TURBULENT VISCOSITY AT NODE POINTS. 
DEFAULT = XMUIMF 

29 

P^RR^Y 

0 

67 

RARRAY(NX) = ANULT, WHERE 

AMULT = COMBINATION OF REMAINING ENTRIES. 

31 

VWALLSTA 

67 

82 

DOWNSTREAM POS. AT WHICH TC 
INJECT TRANSVERSE VELOCITY. 

32 

VWALLVAL 

68 

61 

VALUE CF INJECTED TRANSVERSE VELOCITY. 
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NACA 0015 (Modified) Potentia.] ,F-l£w - 

Data specifications for finite element potential flow analysis over 
single element airfoils has been greatly simplified through use of the COMOC 
grid generator and automated evaluation of the airfoil normal gradient (a3) 
boundary conditions. The basic data requirements therefore, consist of an 
accurate description of the airfoil or viscous boundary shape. The airfoil 
shape is left as a subroutine specification since certain classes such as 
Karman-Trefftz may be generated rather than data specified. Subroutine 
SHAPE is keyed to call various subroutines which generate airfoil shapes. 

An example of data specification for a modified NACA:^0Q15 airfoil is' shown 
in Fig. 2. The airfoil is point specified and the subroutine is called from 
SHAPE when JSHP is set equal to 2 in Namelist NAMEOl. Other airfoil shapes 
may be similarly input by substituting the Fortran subroutine AFSHP in Figure 
2 and specifying JSHP = 2 in NAMEOl. The Fortran deck is placed behind MAIN 
in section, one of the da±a_deck^ — if 

ma y be gene rated from gene ral parameters such as maximum thickness ratio, 
chord length, etc, such as the Karmann-Trefftz class of airfoils, shape 
changes are effected by merely changing the shape generating parameters. 

The required parameters for the Karmann-Trefftz class are trailing edge angle, 
maximum thickness ratio (t/c), angle of attack (a) and camber angle (e). 

These data are input in the Namelist NAME02 and maximum thickness ratio and 
angle of attack are always required since they are used to scale the solution 
domain. Given an airfoil shape, thickness ratio and angle of attack, therefore, 
the appropriate finite element grid for an external flow domain is automatically 
generated and the gradient (a^) boundary conditions are evaluated and applied 
at the appropriate finite element boundaries. 

The following list presents the major sections of the potential flow data 
deck and a brief description of the command and data cards required. 


Section I Fortran Main and Subroutines 

The function of Main is to allocate space for the data arrays. The amount 
of space required is problem dependent and the size of the data arrays must be 
initially guessed at for each different discretization refinement. The actual 
size of the "IZ" array utilized by the program is printed (lARRAY (100) = TREND) 
as illustrated in the output section and the dimensions in MAIN can be reduced 
accordingly on subsequent runs. A list of MAIN for the potential flow test 
case is rllustrated in Figure 2. 

Airfoil shapes are input by subroutine AFSHP which directly follows MAIN 
in the data deck. The maximum airfoil thickness (t/C) must be specified in 
namelist NAME02 to provide the discretizer with the proper domain size scale 
factors. Data requirements consist of the airfoil coordinates along each 
surface non-dimensional i zed by chord length and in the order upper, lower 
proceeding from the leading to trailing edge, and the number of points along 
one surface. 
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c 

-C-G-M-O-C--- 

CCMMON / ARRAYS / 12(90000) 

COMMON / VARBLF / I ARRAY ( 00500 ) , RAP. RAY ( 00500 ) 

— - 

EQUIVALENCE 
CALL ERR SET 

( IARRAY(00092), IZSIZE. ) 

( 207. 500. -1, 1. 0. 217 ) 

1 

00 CONTINUE 


c 

CALL ZEROTK 
CALL RESET 

( 00500. (ARRAY. 0 ) 


CALL RESET 

( 00500, RARRAY, 0.0 ) 

I7STZE = 90000 


CALL RESET 

( TZSTZE, IZ. 0 ) 


CALL BDINPT 
GO TO 100 
END 



2a MAIN Program 


c 

c -- 
c 

c — 
c 

■ ■ c 


c 


i 0‘J 


600 
6 10 


S'uaROuTi ?> ie: afshpid.xiap ,x2Ap,mpts) — — 

PUSPCS6 — PniNf SPEC IF ICATIOM OF, Ah AIRFOIL SHAPE 

CariROIMATES IN CHORD LENGTHS, UPPER SURFACE FOLLOWED 

BY LOWER SURFACE, LEADING EDGE TO TRAILING EDGE, 

100 TOTAL POINTS (HAX.I 

SANE Ni.)r'OER OF POINTS UPPE;R AND LOWER 
MAXIMIJN THICK. RATIO IT/C) MUST G£ SPEC. IN MAMEOl 
DI.'lENSin?.! 0(1), XlAFtl), X2AFII) ' ' ' " 

01 MEN! SI ON X10015I050), X20015I050) 

DATA XI 00 I 5/ . ■ 

* 0., .0023, .0057, .OllA, .0173, .0229, . 03 A3 , ' . 0',5 8 , 

- . 0572 , . 06;>7, .0801 , . 103 , . 1259 , . 1988, .1717, .206, 

N .2518, .3, ..5933, .3591 , .9399, . 9807, . 5265, .5722, 

* ,6295, .6367, ,7939, .8011, .3539, .9156, .9619, 1 . ', 

N ISvO./ 

DATA X20015/ 

0., .0195, .022, .0287, .0336, .0376, .0933, .0973 , 

* ,05. .0532, .05,5, .059, .0632, .0669, .0693, .0721, 

* .0796, .075, .0/99, .0729, .0709, ,0672, .0633, .059, 

Ci .0533, . 0976, .0908, .0335, .0252 , .016, .0077 , '.0, 

7= 18*0./ . . - . - 

NETS = 32 ■ - 

00 100 1=1, NETS 

X lAF ( I ) = X10015 ( I ) 

X2AF ( 1 ) = X2G015I I ) 

WRITE 1 6, 600) ( Xi AF( I ) , I = 1,WPT5) 

V-/RITF (6,610) (X2AF( I ) , 1 = 1 ,.’iPT5) ' " 

RETURN ■ ■ 

FORMAT (//, lOX ,37HOATA SPECIFIED AIRFGIL X 1-C OORD INA T E S , -- 

* I / ,97X, 5E12 .5 ) ) 

FORMAT!/ ,33X,19HX2C00P.0INATES, { / , 97 X , 5E 12 . 5 ) ) 

END ' 


2b Subroutine AFSHP 


Figure 2 MAIN Program and Subroutine AFSHP 
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Section II NAMELIST Specified Control 
and Reference Parameters 


Command 

Name 

Code 

Function 

2DPF 



Initiates COMOC execution 

RENAME 

&NAME01 



Read NAMELIST 6- Set default values 
Fortran NAMELIST integer data 


JSHP 

1 

Karman-Trefftz class airfoil 



2 

Calls user supplied subroutine 
AFSHP to obtain airfoil shape 


NAFPT 


Number of data points specified 
along one airfoil surface. 


NODE 


Greater than or equal to the 
number of solution nodes expected. 


NVRHS 

5 

Dependent variable to be used in 
STRF 


NVP 

5 

Variable vector extracted from 
dependent variable array. 


NIZS 

250 

Data starting point in the IZ 
array 


KDUMP 

1 

Prints data reflection and IZ 
array entry points. 


NMBOUT 

1 

Number of output variables 


NC 


Number of digits to right of 
decimal in print + 3 


LCOL 


Maximum number of columns in 
discretization 


KROW 


Maximum number of rows in 
discretization 


NMOUT 

2 

Prints dependent variable (}> 
solution in node number ascending 
order. 
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Section II NAMELIST (Cont'd) 


2DPF 


Command 

Name 

Code 

Function 

NMOUT 


3 

Geometric form print 

&NAME02 

THKAF 


Airfoil maximum thickness 


ALPHA 


Angle of attack 


BETA 


Karman Trefftz camber angle 


TEANG 


Trailing edge angle 


RNB 


Flow Reynolds number 


COMPX 


Geometric form print compression 
factor (rows). Required if 
NMOUT = 3. 


COMPY 


Geometric form print compression , 
factor (columns) 

ICOND 


. 

Prints NAMELIST data stored in 
arrays lARRAY, RARRAY. 

FEDIMN 



Dimension, Arrays 


Section III Finite Element Discretization 2DPF 


Command 


Functi on 


NX 

NY 

LINKl 9 
LINK2 . 14 


Specifies grid refinement normal 
to the airfoil surface for each 
super element 

Specifies grid refinement along 
airfoil surface for each super 
el ement 

Generates airfoil flow domain dis- 
cretization sized and scaled for 
airfoil in subroutine AFSiRi^ 

Generate vectors for output 
control 
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Section IV Output Specification 


Sample output is listed in the next section 


Command 

COMTITLE 

DONE 

DESCRIPT 204 
DONE 


Function 

Reads title which is printed below the 
COMOC symbol 

Terminates literal data 

Solution print heading 

Same as above 


The next three commands are inter-related and should be fully understood 
prior to changing them. See subroutine BDINPT under subroutines and variables. 


Command 

Function 

DESCRIPT 332 

Solution print parameter titles 
Starting location at RZ(L(32)) 

MPARA -1 

RARRAY locations of multiplier 
to be applied to parameter print 

lONUMB -1 

RARRAY locations of parameters 
to be printed 

DESCRIPT 203 

Titles to head dependent variable 
(<J>) print 

lOSAVE -1 

List of dependent variables to be 
printed. 

lOMULT -1 

List of locations in RARRAY of 
multipliers to be applied to each 
dependent variable for print 


Section V Dependent Variable Boundary and 
Initial Condition Specifications 

The potential flow solution is strictly a boundary value problem, there- 
fore, no initial conditions are required. Gradient boundary conditions are 
internally computed from the specified geometry, and internally applied as in 
equation 16, hence, no boundary condition specification is required. 
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Conmand 


Function 


IPINT -1 Numbers- of-'thei^Dependfeiit and Parameter 

Variable^' in 'the Solution 

Section VI Solution Procedure 2DPF 

Function 

Non-dimensional ize coordinates 
Finite element matrices 
Print node numbers 
Potential Flow Solution (7^(j> = 0) 
Print potential field 

COMOC Print COMOC symbol and title 

specified under command COMTITLE 

EXIT End of execution 


Command 
Link 3 4 

LINK 1 3 

LINK 2 5 

LINK 2 7 

LINK 2 6 


Bradshaw Boundary Layer 

The finite element discretization for 2D boundary layer solutions is one 
dimensional and essentially consists of a column of one dimensional finite 
elements normal to the airfoil surface. Explicit integration of the boundary 
layer equations is performed over the column of elements as it marches in the 
major flow direction. Specification of the finite element geometry has been 
greatly simplified thru the use of geometric progression super elements where 
each element size is a geometric function of the one preceeding it accord- 
ing to the equation 

M+1 

"in < "o * ® j?2 (21) 

In equation (21) p is the specified geometric progression ratio and M is 
the number of finite elements to be generated and scaled by S. For airfoil 
flows, therefore, two super elements are required. The progression for the 
lower surface should be less than unity, while for the upper surface, greater 
than unity to provide discretization refinement in the vicinity of the airfoil 


The following list presents a description of the Bradshaw data d^ck. 
Parameters previously described are"T)m1tted^O“avoM'n^ediMdBTicyT~-ft compTeTe 
listi ng jQLf_jthe_actua-l--data--cleck-4s -in- 4ppendix-B^ — 
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Section II Namelist Specified Control 2DBL 




And 

Reference Parameters 

Command 

Name 

Code 

Function 

FENAME 




&NAME01 

IPTSPL 

1 

T(wall), Patanker & Spaulding 



0 

t(wall), Ludwig Tillman 


ITWALL 

1 

T(wall), du/dy wall 



0 

Use one of the above 


NEQKNN 


Number of variables being integrated 
(must hOit.be greater than NEQ) 


NM 

2 

One dimensional Finite Elements 



3 

Two-dimensional Finite Elements 


NPRNT 


Print page size, in columns 


NIKS 


Number of integral parameters printed 


NTPRNT 

99999 

Suppress integral parameter print 


KNTPAS 


Maximum number of integration steps 
between prints 

&NAME02 

RNULOC 

0. 

Laminar viscosity for Van Driest 
coef. in DFCFBL 



1. 

Turbulent viscosity for Van Driest 
coef. in DFCFBL 


REFL 


Reference length applied to output 


UINF 


Mean Flow Velocity (ft/sec.) 


TOFINF 


Reference temperature (Rankine) 


PINF 


Reference pressure (PSF) 


TO 


Initial station for explicit 
integration in direction of flow (ft.) 


TD 


Distance to final integration 
station (ft.) 


DELP 


Print interval {% of TD) 
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Bradshaw 


Section III Finite Element Discretization 


2DBL 


Command Functi on 

LINK2 14 Forms discretization 

VX2SCL 

ZO, NINTl, Zl, ZPRl, NINT2, Z2, ZPR2 Sample discretization card. 

ZO - first Z coordinate 
NINTl, 2 number of finite 
elements, ZPR1,2 - geometric 
progression ratio 

NDECRD 

NODZ, NOZ, 1, 1, 0, T Sample node selection card 

NODZ first node in discreti- 
zation NOZ number of nodes 
in Z direction 

ELEM Finite element connection 

table 

DONE Ends finite element genera- 

tion sequence 

Section V Dependent Variable, Boundary and 
Initial Condition Specifications 

In this section, initial values of all dependent variables being integra- 
ted must be specified at each solution node. An exception to this occurs when 
integrating the TKE and DISS equations, since a program option exists which 
allows these variables to be internally initialized from an MLT turbulence 
model. Boundary conditions (eqn. 16) are specified in this section in two 
forms, together with a vector of nodes along which each is to be applied. 
Dependent variable values are held constant by simply specifying the dependent 
variable number on the KBNO command card and listing the node numbers in free 
format on cards following, (see subroutine GETBND in the next section). 
Gradient boundary conditions are applied in a special format as described 
under subroutine GETBCD . Other solution parameters may also be input in this 
section in the form of tables to be interpolated as solution progresses (i.e., 
pressure coefficient for dp/dx^ evaluation). Table 1 lists and describes the 
allowed command names and the function of. each. 
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Bradshaw 


2DBL 


Command 
KBNO 1 

BOTTOM 

VXSST 

VPVSX 

VYY 

VYYEND 1 


Command 
LINK2 4 

LINK2 20 
LINKCALL -1 


Section V Dependent Variable, Boundary 

And Initial Condition Specification 


Function 


Nodes where values of the variable 1 
are to be f i xed 

Special command card which fixes all 
nodes along one surface. Other 
options are'Top, Right and Left. 

Tabular input of x, locations; for 
pressure table (ft.) 

Tabular input of edge velocity 

1)2 velocity initialization 

(Non-dimensional input) may be input 
dimensionally if non-dimensionaliz- 
ing factor is specified on the VYY 
card beginning in column 11. A 
convenient method of nondimension- 
alization is to specify the integer 
location of UINF in the RARRAY array. 


Section VI Solution Procedure 


2DBL 


Function 


Calls the continuity equation 
solver 

Bradshaw data initialization 

Subroutines called during the 
integration phase 

1 5 Obtain P&P^ from table 

2 3 Evaluate gradients at airfoil 

boundary 
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5 

6 

Evaluate effective viscosity 


2 

15 

Evaluate the integral parameters 

QKNINT 



Initiates explicit integration 
phase to march the solution from 
TO to TO+TD. At each integration 
station, cal Is are placed to the 
LINKCALL subroutines specified 
under the command (Linkcall). 


Joukowski Wake 

A viscous turbulent solution was performed in the wake of a 12% thick 
JOukowski airfoil at 6° angle of attack. The parabolic Navier Stokes option 
(2DPNS) was utilized and initial velocity profiles at the trailing' edge<were 
experimentally derived (ref. 4). The solution was initialized .01 chord 
lengths upstream of the trailing edge and marched downstream and into the 
wake region. A method of incorporating problem specification changes in 
COMOC such as the sudden ending of the airfoil at the trailing edge is to 
issue a RESTART command which is activated following return from QKNINT. 

Upon restart program variables may be respecified and integration continued. 
For the Joukowski airfoil test case the solution was begun at .'99 chord 
using the 2DBL option with TD set at .01. The restart deck following the 
QKNINT command contained boundary condition information which effectively 
removed the surface gradient. U 2 boundary conditions were applied and LINKCALL 
was modified to perform a 2DPNS solution. The data deck sequence required 
to perform this sequence is described below. Again, duplication of previously 
defined parameters is avoided. A complete listing of the data deck required 
to perform this solution is given in Appendix C. 


Section II Namelist Specified Control and Reference Parameters 2DBL 


Command 

Name 

Code 

Function 

3DBR 



Initiates COMOC executions 

FENAME 



Reads namelist data 

& NAMEOl 

NEQ 


Number of variable, vectors in the 
solution 


LG 


> Number of values in array. CNTPTf.. 


NU2P0S 


Number of values in tables VU2P0S 
And VU2VAL 



Conmand 


Name 


Code 


Function 


2DBL 


NU3P0S 

NEQADD 

NC 

NBC 


&NAME02 

VSTART 


XMUINF 

HSINIT 

C4EDSW 


PRTKE 

PRDIS 

CITKE 

C2TKE 

CKTKE 

CD 

YLTKE 

ESCF 

YPLUS 


Number of values in tables VU3P0S 
And VU3VAL 

Number of variables in IPINT not in 
tegrated until VSTART or C4EDSW 
values of Xj. are reached. 

Number of characters in a print 
word 

Number of a^ , aj boundary 
conditions to be applied. 


% of TD where U 2 is to be 
i ni+fated- when-os^iftg -e OA fT gS - 


Mean viscosity 

Initial integration step size 

location where TKE and dissipa- 
tion become dependent variables 
and begin integrating (see NEQADD 
in NAMEOl) 

TKE from Prandtl mixing length 

Dissipation from Prandtl mixing 
length 

Dissipation production coefficient 
Dissipation dissipation coefficient 
TKE production coefficient 
Constant in dissipation equation. 

Constant for Diss. length calcula- 
tion: foy MET. 

DiSs, with scale factor from MLT 
Non-dimensional Y wall 
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Joukowski Section V Dependent Variable, Boundary 2DBL 

And Initial Condition Specifications 

Command Index Function 

VU3P0S Tabular specification of xi coordinate 

describing geometric boundary variation 
in the Z direction. 

VU3VAL Values of Z at each VU3P0S xi .coordinate 

The 2DBL and 2DPNS equations are scaled 
to allow for grid growth as the solution 
proceeds 

KBNO I Nodes where variable is to be fixed 

1=1 Variable No. 1 - Uj velocity 

I = 2 Variable No. 2 - U 2 velocity 

I = 5 Variable No. 5 - Turbulent Kinetic 

Energy 

1=6 Variable No. 6 - Dissipation Function 

VYY 

VYYEND 1 Only Ui need be initially specified. For 

boundary layers U 2 is computed from the 
continuity equation. The continuity 
equation is initiated by setting VSTART 
in NAME02, and U 2 is counted as a var- 
iable in NEQ but not in NEQKNN in NAMEOl 

TKE and DISS, variables 5 and 6 
respectively, are initialized from the 
MLT model. Setting NE1E2 equal to 2 and 
E1E2SW to a small increment of TO will 
initiate a laminar-turbulent (MLT) 
sequence and NEQADD = -2 will maintain 
TKE + DISS as solution parameters to be 
evaluated from MLT results. C4ED SW in 
NAME02 is the location where varia- 
bles 5-6 become^solution dependent 
variables utilizing the values stored 
during their parameter status as 
initial conditions. 
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Joukowski 


Section VI Solution Procedure 


2DBL 


Command 
LINK2 21 
SAVETAPE 

RESTART 


The following explains the restart deck in Appendix C as used 
specifically for the Joukowski test case. 

Section II Namelist Specified Control 2DPNS 

And Reference Parameters 

Name Code Function ' 

Reads Namelist data 

KDUMP 0 No print of input 

1 Print input as encountered. 

NEQADD 0 All equations integrated 

immediately upon restart 

NE1E2 0 Do not evaluate MLT anywhere 


Command 

NAMELIST 
& NAMEOl 


Index Function 

Obtain pressure from the Cp table. 

10 2 Store the entire data arrays set 

on unit specified at each print 
station. The second integer indi- 
cates the number of back-up writes 
to retain for restart. 

01 To initiate a restart and allow 

changes in the data deck. Change 
equations to 2DPNS 
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Joukowski Section V Dependent Variable 2DPNS 

Boundary and Initial Condition 




Specification 

Command 

Index 

Function 


IBORD Set up vector of boundary-nodes for 

boundary condition input. 


RIGHT 




Freestream node along upper 
surfaces 

DONE 




Ends literal data 

KBNO 


1 


All nodes in solution for Variable 
1 (ui) 

DONE 




Terminates Literal Data. 

KBNO 


2 

1 

Dep. Var. 2 boundary conditions 

BOTTOM 

0,0,0, 

0.0, 

2, 4.321, 2 

0.0,2 sets 41= 0.0 

4.321,2 sets '^'-37 

Same as for bottom. 

TOP 

0,0,0, 

0.0 

2, -4.321, 2 


These cards set the boundary 
conditions required to integrate 
Uo to form 2DPNS equation system. 
±4.321 is a code indicating that 
the as frees tream boundary con- 
dition is to be obtained from 
du/dx. 


DONE 



Ends literal data 

KBNO 

5 


Removes fixed restrictions on 
TKE new airfoil surface 

DONE 



Terminates Literal Data 

KBNO 

6 


Removes fixed restriction on 
dissipation near airfoil surface 

DONE 




I ARRAY 

145 

0 

imi:n = 0 

RARRAY 

15 

7 -1 

H = (non-dimensional) 

RARRAY 

22 

1.001 -1 

Reset final station 


(non-dimensional ) 
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Section VI Solution Procedure (Cont'd) 


2DPNS 


Joukowski 

Command - Index 

LINKCALL -1 T 

1 5 

2 15 

5 6 

3DPNS 

SAVETAPE 10 2 

QKNINT 

EXIT 


Function 

Calls routines to solve the 
2DPNS equations 

Evaluate dp/dx 

Evaluate integral parameters 

Evaluate diffusion coefficients 

Flags 2DPNS solution 

Saves data on unit 10 modulo 2. 

Initiates 2DPNS equation 
integration. 

End of job 


Selected Print Samples 

The output package in COMOC is quite versatile and allows the user 
substantial print format control. The standard print may be described in 
the three distinct classes of problem identifying titles, data reflection 
print and solution print. Each of these operates under control of the user 
as described in the previous section. Three different titles may be input 
under different command names for print below the COMOC symbol on command 
COMOC, print at the beginning of the data reflection or as a descriptive 
heading at the beginning of each solution print. 

Data reflection print is accessed by setting KDUMP = 1 in NAMEOl and the 
print format is illustrated in Figure 3. Each card in the deck is printed 
as it is encountered. In addition, vectors filled by array data are printed 
directly following the data reflection, thus providing a positive verifica- 
tion of proper input. Common data errors such as not specifying "T" or 
blank card delimiters are easily detected since the reader will continue to 
read cards until it encounters one of these. The print subsequently appears 
as continuous data reflection of all the cards under the original command 
name, thus providing a quick and positive data check. Scalar data is stored 
in the arrays RARRAY and I ARRAY. Print of these data is accessed by inserting 
the comnand ICOND. Figure 4 illustrates the print format. This prov^ictes 
a quick check of the namelist data specification and program computed scalars 
which are stored in these arrays. 
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Figure 4a Scalar Integer Data Contained in lARRAY and Passed 
Through COMMON/VARBLE/ , (First 200 Locations) 
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Figure 4b Scalar Real Data Contained in RARRAY and Passed 
Through COMMON/VARBLE/ , (First 200 Locations) 



Solution print may be subdivided into two parts for purposes of discussion. 

The first contains header information which includes narrative titles speci- 
fying the option being executed. This is followed by a titled table of 
parameters as specified in Section IV of the data deck. Figure 5 illustrates 
one of the more elaborate tables which includes mean values of various flow 
and thermodynamic parameters in four different sets of units. Integration 
parameters such as Xj station, current integration step size, cumulative 
number of passes through the derivative evaluation routine, and number of 
print stations are also given to provide information on integration status. 

The header section is followed by print of the solution variables at 
their current computed values. The variables specified under the lOSAVE 
and lOMULT commands in Section IV of the data are printed according to the 
digit format specified in NAMEOI parameter NC. Titles for each variable 
printed are specifiable under command DESCRIPT-203 and appear as illustrated 
in Figure 6. All dependent variables and their derivatives and node specified 
solution parameters may be printed in this form. 

Geometric form print of solution variables is the default option for print 
at each solution print station. This form of print is illustrated in Figure 
6 where the airfoil surface is envisioned at the center of the geometry and 
the discretization proceeds above and below the surface. A more graphic 
example of the usefulness of this option occurs in the potential flow solution 
where a two-dimensional finite element grid is required. Figure 7 illustrates 
that by geometrically ordering the print, values of potential function at 
points of interest with respect to the airfoil are readily identified. In 
order to keep this print form within a few pages, however, compression factors 
must be applied to the data. These are specified as COMPX and COMPY in NAME02, 
and higher numbers provide the most compression. Unfortunately, compression 
causes some of the data to be eliminated from the print. The region of maxi- 
mum loss is in the highly refined and more interesting portion of the flow- 
field. This problem can be overcome, however, by resorting to the more stan- 
dard columnar form of print which is accessed by setting NMOUT = 2 in NAMEOI. 

In addition to standard print of discrete variables at selected print 
intervals, the option exists to print integral solution parameters in the 
2DBL and 2DPNS COMOC options at each integration station. This print, as 
illustrated in Figure 8, provides an indication of solution progress between 
print stations. Solution of these parameters must be requested (2-15) under 
the LINKCALL command in Section VI of the data deck. 

Upon call to QKNINT, a series of prints is initiated. The first, illus- 
trated in Figure 9, prints the variable numbers and types in the solution. 

If a restart unit is requested, the unit number is printed. The order of 
calls for solution process is listed and identified, followed by the vari- 
ables and multipliers to appear following the header page at each DEEP. 
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Figure 5 Typical Solution Header Print 
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Figure 6 Typical Dependent Variable and Derivative Print 
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Figure 7 Node Map of Potential Flow Solution Domain 
NACA-0015 (Mod.), First Quarter Chord 
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Figure 8 Integral Parameter Print 



solution print interval. Following this print, a call to ICOND is initia.ted 
to print the RARRAY and lARRAY scalar lists. This is followed by a node map 
print which matches node numbers with their coordinates. A subsequent 
standard solution print, as described above, is output which lists the 
specified initial conditions. 


4 variables being INTFGKATED, 
15 6 2 

A VARIABLES IN SOLUTION. 
15 6 2 

SAVE r-ITRUT ON UNIT 10 ‘ 

ORDER OF CALLS AT END OF OKNUIN 

LINKK 5) - DPOX 

L1NK2I 4) - CONTES 

LINK2(15) - TRBTHK 
LIMK2I 3) - ViLFLXS 

LINKS! 6)~- SETDIF' 

PRINTOUT VARIABLES 

1248 2248 5248 6243 1247 

1249 . 22^9 52^9 CZ/tJ) 305 ’ ’ 7. 

PRINTOUT VARIABLE MULTIPLIERS. 


Figure 9 Solution Sequence Print, 

Obtained Upon Entering QKNINT 




COMOC COMPUTER PROGRAM SYSTEM 


General Overview 

The COMOC computer program system is rapidly developing into its intended 
design state as a general purpose differential equation solver. Present 
capabilities include various fluid mechanics solution options (2DBL, 3DBR, 

2DNS, 2DPNS, and 2DPF) including various turbulence modeling and thermodynamics 
evaluation. The most recent capability combines (2DPF, 2DBL and 2DPNS) into a 
coupled weak interaction aerodynamic iterative flow solution sequence. The 
program I/O is highly developed and approaching a state of data specified 
mathematical operations. The finite element method is utilized as the basic 
numerical algorithm, thus providing the ability to easily model general geo- 
metric boundaries and allows for simplified boundary condition specification. 

The program is written entirely in Fortinan and presently consists of 
approximately 19000 cards. Array storage is dynamically allocated according to 
problem size, making problem size limitation strictly dependent upon the 
computer size on which it is run. This section gives a brief description of 
the program flow followed by a list defining the subroutines and some of the 
more important variables in the program. Finally, some examples of diagnostic 
print obrainable by code in Namelist NAMEOl are illustrated. 

The COMOC macro-structure is illustrated in the flow diagram of figilre 
10. Main allocates core for the data arrays. The input module consists of the 
sub-routines used to control the program flow and read unformatted data. All 
data is read by subroutine REDREL controlled by READER . Command data are 
interpreted by BDINPT which is the program director. Interpretation includes 
the three categories of program control, integer data and real data. Command 
data, therefore, controls the program sequence and directs the filling of 
arrays. The LINK name is utilized to perform a sequence of operations which 
require several subroutines. This manipulation helps insure that certain in- 
ternal operations are performed in the proper sequence at the expense of some 
user control. Vector initialization and discretization are user controlled, 
but must be performed prior to the integration phase. Integration is initiated 
by calling subroutine QKNINT and integration progresses via the integrator 
QKNUIN . Derivative vectors to be integrated are formed in DERVBL . The user 
has control over the equations evaluated through specification of dependent 
variables in the IPINT array and NEQKNN in NAMEOl. User control of the sub- 
routines to be invoked each integration step is exercised under the LINKCALL 
command. Output is controlled by sub-routines REOUTP (prints node map) and 
FEOUTP (prints data array) which are accessed through calls to LINK2-5 and 
LINK2-6 respectively. Execution proceeds until T = TF whence control is 
returned BDINPT . Command END is used to terminate a data case and EXIT 
terminates execution. 

The following pages list and give a brief description of the subroutines 
and some of the more pertinent variables in the program. 
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Figure 10 COMOC Macro-Structure 










SUBROUTINE DESCRIPTIONS 


TEE FOLLOWING PAGES CONTAIN A BRIEF DESCRIPTION OF THE SUBROUTINES 
IN TEE CCRCC CC^'fUTEF PROGRAM. 


NAEES IN PAPENTEESF.S INCICATE CALLING ROUTINES. 

IF NO NAME IS ENTERED THEM SEVERAL ROUTINES PLACE CALL. 

MAIN 

THIS IS TEE MAIN CCNTRCL PROGRAM WHICH INITIALIZES THE RARRAY, 
lAPPAY AND tee IZ ARRAYS TO ZERO. 

TO CEANGE THE VARIABLE SOTRAGE CAPACITY CF THE IZ ARRAY, RESET 
THE CIEEESICN CF 'll' ANC, ACCORDINGLY, THE VALUE OF 'IZSIZE'. 

AFTER INITIALIZATION TEE CCNTRQL ROLT INE BDINPT IS CALLED. 

ECINPT (MAIN) 

THIS IS TEE CCNTRCL FCLTINE WHICH INITIALIZES VECTORS AND CONTROLS 
THE FLOW OF THE PROGRAM ACCORDING TO LSER INPUT. 

( SEE TABLE 1 FOR INPUT DESCPIFTICN. ) 

LINKl 

PLACE CALLS TO TEE FCLLOkIN'G ROUTINES. 

2. NCEcLM 

3. GECMFL 

A. GETPPR, GT'JEDG 
5. GETPPR, FRSGRC 

7. eCCMDT 

8. CERVBL 

L INK2 

PLACE CALLS TC TEE FCLLCWING ROUTINES. 

1. OFCFNS 

2. CFCFBL 

3. WLFLXS 
A . CONTES 
5. PECUTF 

•6. FECUTP 

IF ' ITEM' .GT. 0, WRITE 'PLOTS' DATA ON UNIT 'ITDB'. 

IF RESTART CODE 'NRTAFE' IS . GT . C, WRITE RESTART DATA ON TAPE 'NRTAPE'. 
7 . STRF 

S. CRFCQL, CRECNS, CRFOGS 

10. T BLIMP 

11. CFFPEI 

12. H2M1X 
13 . XYCPDM 
lA. DSCRTZ 
15. TRETHK 

21. PFRMCF 

22. ARFOIL 

23. SYMELM 
25. BLTINT 
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LINK 3 


FL/iCE U>LLS TC TFF FCLLCWIKG ROUTINES. 

1. NUNORV 
2 . RITE 

4. OPCXTe, CI^'E^ .4^jC STORE DEP . VfiR. LOCATION LIST IN IZ(INPINT) 

5. LTH 

L INKA 

PLACE CALLS TO TEE ECLLCwING ROUTINES. 

2 . QK^U!^ 

9. PCTENT 

LINKS 

PLACE CALLS TO TEE FCLLOV^ING ROUTINES. 

1. NWGECE 

2. CALL RCUTINES FROP LINKCALL LIST AT END OF QKNUIN. 

3. CPINIT 

5. SCHPRN 

6 . S ETOI F 


AESAVE 

CCMFUTE TEE SUN OF AESCLOTE VALUES OF A SEQUENCE OF NLINSERS. 

AODDEL ( ELFM, GETRCO, GETONO, SETUP ) 

ACO CP DELETE ENTRIES IN AM INTEGER ARRAV DEPENDING ON THE 
VALUE CF ' KT VP' . 

KTYP = 1, CELETE 
KTYF = 2, ADC 
ASHSQ (STRF) 

eCCLEAN ASSEPCLY OF SQLARE SYNMETRIC MATRIX. 

ASMVeC 

EOOLEAN ASSEMELY OF COLUMN MATRIX. 

AVRG 

COMPUTE TEE arithmetic AVERAGE OF 'NUPO' ENTRIES IN AN ARRAY. 

EANCFG (STRF) 

SYMPETPIC EANCEC CFCLESKY LINEAR ALGEERAIC E13UAT1GM SOLVER. 

8CGNDT { LINKK7) ) 

GRADIENT RCUNCiRY CCNDITICN SPEC I F I C AT IC N (A1,A3) FOR DEP. VARIABLE. 

-1.23A = USE DUDY FOR A3. 

+ OR -A. 321 = SET DVDY = ~ DUDX 

BNOSET ( GETGCO .GETBNO I 

DETERMINE NOCES TO PE INSERTED INTO BOUNDARY ARRAY. 

CALORD (eOINPT) 

PRINT CALL CRDER OF PGUTINES USED FOR VARIABLE PARAMETERS. 

PRINT INTEGRATIGN VAPIAELE NCS . ANC CEPENOEMT VARIABLE NOS. 

PRINT LIST OF PARAMETERS TO BE PRINTED IN OUTPUT ROUTINE. 

CELS (CSCPTZ) 

CCMPLTE THE NU'*CER CF COLUMNS, 'LCCL', IN THE OUTPUT DISPLAY 
AND SET UP TEE FOLLCWING ARRAY, 

INCCLU) - NC. CF NCCES IN COLUMN J. 

CCMCC (ECIMPT) 

PRINT THE CCMCC SYMBCL CN TEJO PAGES ALONG WITH ASSOC. TITLE CAROS. 
CONTES IL1NK2IA)) 

RUNNING SMOOTH CONTINUITY EQUATION SOLVER TO COMPUTE U2 UP COLUMNS 
CF NCDES AFTER VSTART EAS BEEN REACHEC. 

CPINIT IDIMEN) 

CCMPUTE CP IMF AT TS INF . 
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DELADO CAODDEL) 

ADD ENTRIES TO AN INTEGER ARRAY *NSIOE» AT A TIME. 

OELELM (DELNOD) 

DELETE ENTRIES IN AN INTEGER ARRAY 'NSIDE* AT A TIME. 

DELETE IDSCRTZJ 

DELETE NODES THAT ARE NOT CONNECTED TO ANY ELEMENTS. 

DELNOD (ADDDEL) 

SET UP CALL TO OELELM AND SUPPRESS ZERO ENTRIES IN A VECTOR. 

DERVBL ILINKlISn 

fORM THE DERIVATIVE OF THE ORDINARY DIFFERENTIAL EQUATION FIRST 
ON UX-VELOCITY IGLOBAL CONTINUITY) AND OTHER DEPENDENT 
VARIABLES INCLUDING SPECIES CONTINUITY, ENERGY, LONGITUDINAL 
AND LATERAL MOMENTUM, IF REQUIRED. 

DERVDX (CONTES) 

3-POINT FORWARD OIFF. FORMULA TO COMPUTE DERIVATIVE IN DOWNSTREAM DIR. 
DESCRPIKA) IBOINPT) 

DPDXTB 

COMPUTE DPDX FROM PRESSURE TABLE. OPOX IS CONSTANT OVER INTERVALS. 
DFCFBL ILINK2<2)) 

COMPUTE TURBULENT VISCOSITIES FOR DEPENDENT VARIABLES FROM 

1. TKE - DISSIPATION EQUATIONS OR 

2. MIXING LENGTH THEORY, OR 

3. COMBINATION OF BOTH TKE - OISS. AND MLT. 

OFCFNS (LINK2(1J) 

COMPUTE LAMINAR VISCOSITY ACCORDING TO TEMP. AT NODES 
USING SUTHERLAND'S LAW. 

OIMEN ILINKSUn 

COMPUTE NON-OIMENSICNALIZING FACTORS USED IN PROGRAM. 

DRHOBL ILINK2C?)) 

CALL IF IGAS = 0 IN NAMEOl 

COMPUTES THE TEMPERATURE AND DENSITY USING A SIMPLIFIED 
ENERGY EQUATION. ISOENERGETIC MIXING FLOW WITH 2 SPECIES. 

ORHOGS (LINK2(g)) 

CALLED IF IGAS = 1 IN NAMEOl 

COMPUTES THE TEMPERATURE, DENSITY AND SPECIFIC HEAT ON A NODAL 
BASIS AS A FUNCTION OF PRESS., ENTH. , VEL. AND SPECIES COMPOSITION. 
IF NGETH = I IN NAMEOl, THE FIRST PASS THROUGH THIS ROUTINE 
WILL RETURN ENTHALPY WHEN GIVEN THE TOTAL TEMPERATURE AT THE NODES. 
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CSCRTZ (LINKZa^l) 


SE^" UP SGLUTICN CCNAIN C ISCRET IZAT ION . 

SET UP SCALE FACTORS, GENERATE NODES AMO ELEMENTS, COMPUTE ROW 
ANC COLUMN KEYS FOF SCLUTICN AMO PRINTOUT, AND SCALE INPUT 
COORDINATES 6V RGFL. 

TFE FOLLOWING KEYWCRCS ARE USED TO TRANSFER FLOW - 


KEYWORD 

FCRMAT 

ELCCK CCLS. 

CESCRIPTION 

VXISCL 

A8 

L - 8 

R6AC the FOLLOWING IM FREE FORMAT 
XO - START OF XI GEOMETRY. 

NDIVl - NO. OF DIVISIONS IN FIRST I NT . 
XI - XI POS. AT END OF INTERVAL. 

PR I - PROGRESSION RATIO FOP SPACING. 

NOIV2 - NO. OF DIV. IN 2ND INTERVAL. 

X2 - X2 POS. AT END OF 2ND INTERVAL. 

PR2 - PROG. RATIO FOR SPACING. 

. . . CONTINUE WITH ND1V3 X3 PR3, ETC. 
UNTIL A SCAN DELIMETER 'T' OR A BLANK 
CARC IS EMCOUNTERRED. 

VX2SCL 

A8 

1 - 8 

SAME AS VXISCL BUT FOR DIRECTION 2. 

NDSCRO 

A8 

1 - 8 


N 1 

FRE5 

AFTER COL. 8 
FOR N1 E( 

FOR Ml N£ 

)N NOECPO CARD. 

3 -1, A RECTANGULAR GRID IS FORMED 

FRCM THE VXISCL AND VX2SCI. INPUT. 

-1, READ IN RECTANGULAR MESHES, A PER 
CARD IN FREE FORMAT. 

TERMINATE READ WITH A BLANK CARD. 
NRL, NRU, NCL, NCU 
MU - LCWER ROW NUMBER 
MRU - UPPER ROW NUMBER. 

NCL - LOWER COLUMN NUMBER. 

NCU - LPPER COLUMN NUMBER. 

E . G. 

3 8, 1 6 

PRODUCES A MESH OF NODES FROM ROW 3 
THRU ROW 8 AND FROM COLUMNS 1 THRU 6 
WITH SCALE FACTORS GENERATED ACCORDINGLY. 

ELEM 

A3 

1 - 8 

CENERATF ELEMENTS FROM RECTANGULAR MESH. 

NL 

FREE 

5FTER COL. 8 
:n ELEM CARD 

+l = ACC ELEMENTS IN ELEM. 

N2 

FREE 

IFTSR M 1 

-1 = delete ELEMENTS IN ELEM. 

N3 

FREE 

6FTER N2 

TURN DIAGONALS BELOW ROW N3. 

CCNE 

A8 

1 - a 

SCALE XICOR BY XSCALE 

SCALE X2COR BY YSCALE AND RETURN. 


ANY OTHER KEYwGRC ENCCUNTERREO HILL CAUSE A RETURN FROM DSCRTZ. 


« * * mote ^ t VXISCL IS OPTIONAL FOR ID ELEMENTS. 
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DUDY IWLFLXS) 

THREE POINT INTEGRATION FORMULA FOR COMPUTING OUDY. 

ELEM lOSCRTZ) 

GENERATE ELEMENTS AS A FUNCTION OF NODE COORDINATE INPUT. 

USED PRIMARILY FOR A RECTANGUALR DOMAIN. 

ELKEY2 (L7HJ 

GENERATE KEYS CORRELATING ELEMENT DOF TO. SYSTEM DOF. 

FEDIMN (BOINPT) . 

SET UP DIMENSIONS OF VARIABLE LENGTH ARRAYS USED IN THE 
SYSTEM- FINOS LOCATIONS OF OUTPUT ARRAYS FOR 'FEOUTP*. 

IF *KDUHPV = 1, print LOCATION OF ENTRY POINTS IN 'I2* ARRAY. 

FENAME (BOINPTi 

THIS ROUTINE CONTAINS A LIST OF ALL EOUIVALENCED VARIABLES 
IN THE I ARRAY AND RARRAY VECTORS. 

MOST DEFAULT VALUES ARE ALSO SET IN THIS ROUTINE. 

CALL NMELST to READ IN NAMEOl AND NAME02 NAMELISTS. 

FEPLOT (STOUTIJ 

GENERATE DATA TO BE USED FOR PLOT PACKAGES. 

FINOBE IBOINPT) 

DETERMINE A SERIES OF BOUNDARY ELEMENTS AS A FUNCTION OF 
IMPUTING BOUNDARY NODES IN COUNTER-CLOCKWISE ORDER. 

ON FIRST PASS. IF IBORD IS READ, FIND BORDER ELEMENTS 

AND REORDER NODES SO THAT FIRST TWO ARE ON THE BOUNDARY. 

GENDA 

EXTRACT AN ELEMENT VECTOR FROM A GLOBAL VECTOR USING THE ELEMENT 
CONNECTION TABLE ‘INODE'. 

GEOMOR IGSOMFLI 

COMPUTE ENTRIES FOR B112 AND 8113 MATRICES. 

GEOMFL ILlNKlian 

GENERATE THE UNIQUE ELEMENT MATRICES AND VECTORS. 

SET UP THICKNESS VECTOR ITK. 

GENERATE LENGTH « THICKNESS ARRAY IXIP2. 

GENERATE AREA ♦ THICKNESS ARRAY lAREA. 

GENERATE 8112 MATRIX. 

GENERATE 8113 MATRIX. 

IF KOOG .GT. 0, PRINT ELEMENT NO., NODES OF ELEMENT AND 
COORDINATES OF NODES FOR ELEMENTS FROM *IBOT‘ TO MTOP*. 

AFTER THE ELEMENT LOOP IS COMPLETED, PRINT THE VECTORS AND MATRICES 
THAT WERE GENERATED IN THE ELEMENT LOOP. 

GETAOO 1 FEDIMN J 

IBM 360 ASSEMBLER LANGUAGE ROUTINE TO GET MACHINE ADDRESS OF VARIABLE. 
GETALC (XYCRDMJ 

COMPUTE ‘ALC AS THE SHORTEST SIDE OF ALL THE ELEMENTS IF IT IS 
NOT READ IN NAME02. 
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GETBCO IBDINPT) 

INPUT IS SIMILAR TO GETBNO EXCEPT THERE IS ONLY ONE BLOCK PER CARO - 


KEYWORD 

FORMAT 

BLOCK 

COLS. 

DESCRIPTION 

(SAME) 

A3 

1 - 

8 

SAME OEF. AS IN GETBND. 

KEYWORD 

FORMAT 

BLOCK 

COLS . 

DESCRIPTION 

KODEl 


AFTER 

8 

SAME AS KODEl IN GETBNO. 

K00E2 


AFTER 

KODEl 

SAME AS K00E2 IN GETBND. 

K00E3 

FREE 

AFTER 

K00E2 

SAME AS K0DE3 IN GETBND. 

A1 

FREE 

AFTER 

K0DE3 

VALUE OF A1 FOR THIS BOUNDARY. 

MAI 

FREE 

AFTER 

A1 

RARRAY MULT. FOR Al. 

A3 

FREE 

AFTER 

MAI 

VALUE FOR A3 BOUNDARY CONDITION. 

MA3 

FREE 

AFTER 

A3 

RARRAY MULT. FOR A3. 

♦MAI = MULTIPLY BY RARRAY(MAl) 
-MAI = DIVIDE BY RARRAY(MAl) 


GETBCM (GETBCO) 

extract BOUNOARY CONDITION VECTORS FROM INPUT DATA. 

GETBNO (BDINPT) 

establishes the BOUNDARY NODE VECTOR FOR EACH DEP. VAR. USING 
EITHER THE WORD ‘ADD' OR SIMPLE GEOMETRY OF THE PROBLEM WITH THE 
FOLLOWING KEYWORDS AND COOES - 

EACH CARO IS DIVIDED INTO FOUR IDENTICAL BLOCKS OF 20 COLUMNS EACH. 
ALL BLOCKS ARE OF THE SAME FROHAT SO THAT A 
DESCRIPTION OF ONE BLOCK ONLY WILL BE GIVEN. 

THE BLOCKS START IN COL. 1, 21. 41 AND 61. 


KEYWORD 

FORMAT 

BLOCK COLS. 

DESCRIPTION 

TOP 

A8 

1 - 


ACROSS TOP FROM LEFT TO RIGHT. 

-TOP 

A8 

1 - 

8 

ACROSS TOP FROM RIGHT TO LEFT. 

BOTTOM 

A8 

1 - 

8 

ACROSS BOTTOM FROM LEFT TO RIGHT. 

-BOTTOM 

' A8 

1 - 

8 

ACROSS BOTTOM FROM RIGHT TO LEFT. 

RIGHT 

A8 

1 - 

8 

UP RIGHT HAND SIDE. 

-RIGHT 

A8 

1 - 

8 

DOWN RIGHT HAND SIDE. 

LEFT 

A8 

1 - 

8 

UP LEFT HAND SIDE. 

-LEFT 

A8 

1 - 

8 

DOWN LEFT HAND SIDE. 

(BLANK) 

A8 

1 - 

3 

IGNORE BLOCK. 

ADD 

A8 

1 - 

8 

CALL ADODEL TO INSERT ENTRIES. 
IGNORE BLOCK COLS. 9 - 20. 

DELETE 

A8 

1 - 

8 

CALL AODOEL TO DELETE ENTRIES. 
IGNORE BLOCK COLS. 9 - 20. 

DONE 

A8 

1 - 

8 

LEAVE ROUTINE. 


FOR THE FOLLOWING KEYWORDS, THE THREE CODES (WE‘LL CALL THEM 
KODEl, K00E2 AND K00E3 FOR CONVENIENCE) WILL DETERMINE WHICH 
NODES WILL BE SELECTED. 



FORMAT 

BLOCK 

COLS. 

DESCRIPTION 

KODEl 

14 

9 - 

12 

ROW OR COLUMN DI SPLACEMENT FROM 





EDGE BEING DESCRIBED (DEF. = 0). 

K0DE2 

14 

13 - 

16 

POS. IN LINE TO START (DEF. = FIRST). 

K0DE3 

14 

17 - 

20 

POS. IN LINE TO END (DEF. = LAST ). 
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GETFSL 1 DFCFBL.TRBTHK) 

FIND BOUNDARY LAYER THICKNESS, DELTA, AND NODE AT WHICH IT OCCURS. 
GETINO (GEOMFL) 

EXTRACT NODE INDICES FOR ELEMENT BEING PROCESSED. 

GETPPR ILINKIIBI) 

TABLE LOOK-UP OF PRESSURE AND OPOX AS FUNCTION OF DOWNSTREAM STATION. 

GMADD 

GENERAL MATRIX ADDITION. C = A + B 
GPAHFT I THERMO! 

MULTIPLE SPECIES THERMODYNAMICS. 

IF NGETH .GT. 0, COMPUTE ENTHALPY DISTRIBUTION FROM TOTAL TEMP. 

IF NGETH .LE. 0, COMPUTE TEMPERATURE, DENSITY, SPECIFIC HEAT 
AND MACH NO. AS FUNCTION OF PRESSURE, ENTHALPY AND SPECIES COMP. 
H2MIX {LINK2I12)) 

COMPUTE THE MIXING EFFICIENCY HRSOOT AND THE MASS FLOW HOOT. 

ICOND IBDINPT! 

PRINT INTEGER AND REAL INITIAL CONDITIONS. 
lARRAYIl! - IARRAY(4001 
RARRAY(l) - RARRAYI400) 

LOC (MPRD) 

COMPUTE VECTOR SUBSCRIPT FOR AN ELEMENT IN A MATRIX OF 
SPECIFIED STORAGE MODE. 

LOCATE 

FIND THE LOCATION OF *M» IN THE ARRAY *NA» AND STORE IT IN *N*. . 

LOOK 

LINEAR INTERPOLATION ROUTINE. 

LSFT (MISOIVJ 

GENERATE A LEAST SQUARES FIT THRU A SERIES OF POINTS. 

L7H (LINK3I5I) 

COMPUTE ELEMENT LINKING KEYS FOR SOLVER ROUTINE 6ANCH0. 

MATSUM 

COMPUTE Am = BID COEF * CD), i = 1, N 
MINMAX IDSCRTZ, ORDER, SETSCL) 

COMPUTE THE MINIMUM *MN* AND MAXIMUM 'MX* ENTRIES IN AN ARRAY 
AT LOCATIONS MMN* AND 'IMX' IN THE ARRAY. 

MISDIV ICONTES) 

POLYNOMIAL FIT THRU *NPT* POINTS OF THE M-TH ORDER. 

NPT MUST BE AN ODD NUMBER. 

MNMX (EL EM) 

FROM AN INTEGER VECTOR I NA CONTAINING NN ENTRIES, STORE THE FOLLOWING 
LOW - POSITION IN MNA* OF MINIMUM. 

LHI - POSITION IN 'INA' OF MAXIMUM. 

MN - MINIMUM VALUE IN MNA’. 

MX - MAXIMUM VALUE IN ‘INA* 

MPRD 

MULTIPLY TWO MATRICES AND STORE IN RESULTANT MATRIX. 

C = A >t: B 
MTRA I GEOMFL) 

FIND THE TRANSPOSE OF A GENERAL MATRIX. T = A TRANSPOSE. 

NBNDRY ILINKSm) 

THE VALUE OF 'NBScT* DETERMINES THE OPERATION OF THIS ROUTINE, 

NBSET = 1, SET UP INTEGRATION NODES AND STORE DEPENDENT VARIABLE 
INTO 'YY* VECTOR. 

NBSET = 0, SET UP INTEGRATION NODES AND RETRIEVE DEPENDENT VARIABLE 
FROM ‘YY* VECTOR. 
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NDECRD (OSCRTZ) 

GENERATE NODE COORDINATES FROM SUPER ELEMENT INPUT DATA. 

NMELST (FENAME. BDINPTl 

READ NAMELISTS NAMEOl AND NAME02. 

NAMEOl = INTEGER INPUT. 

NAME02 = REAL INPUT. 

NODELM (LINKlian 

SET UP THE ARRAYS 'lELS* AND ‘lELEH* TO STORE THE NUMBER OF 

ELEMENTS PER NODE AND A LIST OF ELEMENTS CONNECTED TO EACH NODE. 
ALSO AVTHK AND AVAREA 
NWGEOM (LINKSmi 

COMPUTE H21. G22, G23 AND FI FOR VARIABLE GEOMETRY PROBLEM. 

ORDER ICOLS.ROWS.XYSCAL) 

ORDER 3 ARRAYS ACCORDING TO THEIR XI AND X2 COORDINATES. 

THE THIRD ARRAY WILL CONTAIN THEIR ARRAY LOCATIONS. 

OUTNOD 

PRINT AN INTEGER ARRAY ALONG WITH A 32 CHARACTER TITLE. 

OUTPG (GEOMFL) 

PRINT THE ELEMENT NO. AND MODE CONNECTIONS AND NODE 
COORDINATES FROM THE GEOMETRY ROUTINE ‘GEOMFL*. 

OUTVEC 

PRINT A REAL ARRAY ALONG WITH A 32 CHARACTER TITLE. 

PBLANK (REOUTP) 

INSERT BLANKS IN THE OUTPUT VECTOR ‘P*. 

PFRMCP (LINK2(2U) 

COMPUTE PRESSURE TABLE AND DPDX TABLE FROM CP INPUT. 

PLILNK (REOUTP) 

CONVERT A FLOATING POINT NUMBER INTO *A* FORMAT. 

POLY IMISOIV) 

FUNCTION TO GENERATE COEFFICIENTS Cd) IN Y = Ctl) * X**M 
PRINT A 

PRINT A LIST OF REAL NUMBERS IN ‘A‘ FORMAT. 

PRATIO IDSCRTZ) 

COMPUTE NODES USING PROGRESSION RATIO AND END POINTS. 

PRSGRD 

COMPUTE AXIAL PRESSURE GRADIENT. 

OKNINT (3DINPT) 

INTEGRATION CONTROL ROUTINE TO TRANSFER CALL TO OUTPUT 
PACKAGE AT PRINT STATIONS. 

OKNUIN ILINKA(2) ) 

INITIALIZE INTEGRATION CONSTANTS DURING FIRST PASS. 

COMPUTE STEP SIZE AND NEW VALUE OF DEPENDENT VARIABLES. 

COMPUTE UPDATED PARAMETERS BASED ON UPDATED VALUES OF DEP. VAR. 
OMCONC 

COMPUTE A ROUGH APPROXIMATION OF THE AREA OF SPECIES 

CONCENTRATION AND THE MASS DEFECT XMSDF = ROUALC * lAREA-XSUM), 
WHERE XSUM = AMOUNT OF SPECIES PRESENT. 

READER 

READ INFORMATION FROM SPECIFIED INPUT UNIT. 


REAOVl 

READ FROM INPUT UNIT ACCORDING TO THE FOLLOWING FORMAT, 
COL. 1 - 8, A8 

COL. 9 AND POLL. FREE FORMAT INTEGER OR REAL. 
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RECIP (STOUTl) 

INVERT INPUT VECTOR. 

REDREL IBOINPT.ADDDELl 

SCAN AN 80 CHARACTER CARD IMAGE AND CONVERT THE INFORMATION THEREON 
INTO REAL OR INTEGER NUMBERS ACCORDING TO THE FOLLOWING FORMAT, 
DELIMETERS ARE BLANKS, COMMAS OR COLUMN 80. 

INTEGERS CONTAIN NO DECIMAL POINTS 
REALS CONTAIN DECIMAL POINT AND MAY BE *E* FORMAT 
IF A 'C* APPEARS ON A CARO, THE FOLLOWING CARD IS CONSIDERED A 
CONTINUATION CARD 

SCAN TERMINATES WITH A 'T* OR A BLANK CARD. 

A VALUE MAY BE REPEATED, SUCH AS 20 ^ 5.0 MEANS A SERIES OF 20 5.0 
AN ARITHMETIC PROGRESSION MAY BE INPUT, 

5*1-2. 0 3.0 MEANS 3.0 l.O -1.0 -3.0 -5.0 

A GROUP OF NUMBERS MAY BE REPEATED 
41 2.5 5 1.02 -1.2E-5 MEANS 

2.5 5 1.02 -1.2E-5, 2.5 5 1.02 -1.2E-5, 

2.5 5 1.02 -1.2E-5, 2.5 5 1.02 -1.2E-5, 


REORDR (FINOBEI 

REORDER THE NODES OF AN ELEMENT SO THAT THE FIRST TWO WILL BE 
BOUNDARY NODE SPECIFICATION IN THE MBORO' VECTOR MUST BE 
COUNTER-CLOCKWISE. 

REOUTP (LINK2I5I) 

PRINT THE ARRAY GEOMETRY AND NODE NUMBERS IN A PATTERN THAT 
RESEMBLES PROBLEM GEOMETRY. 

(FEOUTP) 1LINK2I6JJ FEOUTP IS AN ENTRY POINT IN REOUTP. 

PRINT OUTPUT PARAMETERS IN A PATTERN THAT RESEMBLES PROBLEM GEOMETRY. 

IF MAX. SCALE FACTOR EXCEEDS •NSM* (DEF, = 10), TERMINATE THE PROBLEM. 
IF OUTPUT PRINT NO. 'KOUNT', EXCEEDS PRINT LIMIT 'LPRINT* IDEF. = 100) 
TERMINATE THE PROBLEM. 

RESET 

RESET ‘NN* ENTRIES OF ARRAY 'A* TO THE VALUE 'V*. 

RESTOR 

REDEFINE A OEP. VARIABLE IF SOME ENTRIES ARE CHANGED 
WITHOUT INTEGRATION OR ITERATION. 

RITE ILINK312)) 

COMPUTE * NUMBER = (KEY-1 )*10 + NMB *, 

GO TO STATEMENT ACCORDING TO VALUE OF ‘NUMBER*. 

IF 'NUMBER* IS OUT OF RANGE, WRITE TITLE INFORMATION. 

ROWS (DSCRTZ) 

COMPUTE THE NUMBER OF ROWS, • KROW* , IN THE OUTPUT DISPLAY 
AND SET UP THE FOLLOWING ARRAYS, 

INROWll) - NO. OF NODES IN ROW I, 

INDRWdJ - COLUMN NUMBERS OF NODES IN ROW I. 

INDEX(J) - ROW NUMBERS OF NODES IN COLUMN J. 

NOCOL(I) - STARTING COLUMN NO. FOR ROW I. 

RSTRHS {0ERV8L1 

FOR 3DBR U1 VELOCITY WITH FIXED WALL, ALLOW FOR TRIANGULAR 
NODES WITH ONLY 1 OR 2 ELEMENTS ATTACHED. 

SCALEV (FEOUTP) 

CALL SCALE ROUTINE FOR UP TO 10 OUTPUT VARIABLES AT A TIME. 
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SCHPRN {LrNK5(5)J 

COMPUTE THE SCHMIDT AND PRANOTL NUMBERS ON A NODE BASIS. 

SETDER IDERVBL) 

COMPUTE 0* = RHS ! LHS, I = 1, NNODE 
SETDIF ILINK5I6]) 

COMPUTE EFFECTIVE VISCOSITY FOR DEP. VAR. USING DFCFNS AND DFCFBL. 
SETSCL (SCALEVI 

SET SCALE FACTOR FOR AN ARRAY OF REAL NUMBERS AND NORMALIZE THE ARRAY. 
SETVAL 

COMPUTE Ain = BI I) C + 0 
SORT ILINKSIBM 

SORT A VECTOR IN ASCENOING OR DESCENDING ORDER. 

+NN = SORT IN ASCENOING ORDER. 

-NN = SORT IN DESCENDING ORDER. 

STOUTl IFEOUTP) 

DIMENSIONALIZE OUTPUT VARIABLES FOR DISPLAY PURPOSES. 

STRF ILINK2(7)1 

IMPLICIT EQUATION SOLVER SETUP ROUTINE. 

SUMKEY IL7H) 

SET UP NODE KEYS FOR EQUATION SOLVER BANCHO. 

SUTHLD (OIMENI 

COMPUTE VISCOSITY USING SUTHERLAND'S VISCOSITY LAW FOR AIR. 

TAUW (WLFLXS) 

COMPUTE SKIN FRICTION USING PATANKER AND SPALDING OR LUDWIEG - TILLMAN. 
THERMO (ORHOGSJ 

INITIATE CALL TO GPAHFT. 

TRBTHK <LINK2(15)I 

COMPUTE AND PRINT INTEGRAL PARAMETERS. 

IF »ITOA' .GT. 0, WRITE PLOT TAPE FOR INTEGRAL PARAMETERS. 

VARMAX I FEOIMN.FEPLOT ) 

FOR 4-NN, FIND MAXIMUM VALUE IN VECTOR, 

FOR -NN, FIND MINIMUM VALUE IN VECTOR. 

VECFUL IDERVBL) 

MULTIPLY A FULL HYPERMATRIX BY A VECTOR OF LENGTH NN. 

VECMAT IDERVBL) 

MULTIPLY A XYMMETRIC HYPERMATRIX BY A VECTOR OF LENGTH NN. 

VECTA 

BOOLEAN ASSEMBLY OF AN ELEMENT VECTOR INTO A GLOBAL VECTOR 
USING INTEGRATION NODE SEQUENCE. 

VHOLES ILINK2I2^)) 

TRANSVERSE VELOCITY ANO/OR SPECIES INJECTION THROUGH A POROUS WALL. 
ITRANSIENT BOUNDARY CONDITION). 

WLFLXS ILINK2I3)) 

COMPUTE THE SKIN FRICTION DISTRIBUTION AND HEAT TRANSFER 
DISTRIBUTION ALONG THE WALL. 

XYCRDM ILINK2U3)) 

GENERATE VECTORS FOR GRID OUTPUT ROWS AND COLUMNS, 

SCALE COORDINATES WHEN RUNNING VARIABLE GEOMETRY. 

XYSCAL IDSCRTZ) 

COMPRESS A VECTOR OF NUMBERS 'XL* BY SCALE FACTOR 'SCFT*. 

FIND • XYO = MAXIXI) - MINIXl) « SCFT • 

IF TWO ADJACENT POINTS OF ARRAY *Xl* ARE WITHIN ‘XYD* OF EACH 
OTHER, SET THE UPPER VALUE EQUAL TO THE LOWER VALUE. 
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AIRFOIL GFNERATOP AND INITIALIZER. 


The FOLLOWING SLGRCLTINES ARE USED TC, GENERATE AIRFOIL 

CISCRETIZATICN AND IN I T I AL I Z AT I C N FOR POTENTIAL SOLUTION, 


AFSI-P (SHAPE) 

USER SU0PCUTINF TO CEFIM AIPFCIL COCRCINATES. { JSHP = 1 IN NAMEOl) 

AIRFOIL THICKNESS A.NC ANGLE OF ATTACK RUST BE SPECIFIED IN NAME02. 
APFCIL (SYNELP) 

DEFINES SUQPEGICA CATA FOR POTENTIAL FLOW CALCULATION AND MODIFIES 
GPaciENT RCUNDARY CONCIticNS ALONG THE AIRFOIL SURFACE FOR 
ANGLE CF ATTACK. 

BND (REFINE) 

APPLIES BCUNCAPY CCNCITICNS TO REFINEC GRID BOUNDARIES. 

CARRY (RFF INE ) 

EXTRACTS COOPCINATE ANC PARAMETER CATA FOR A PARTICULAR 
SUBREGION FRCP SUEREGICN DATA. 

CPFPHI (LINK2( 11) ) 

CCPPUTES PRESSURE COEFFICIENT ALONG VISCOUS BOUNDARY FROM A GIVEN 
POTENTIAL SOLUTION. 

ElKEY (REF INE ) 

generates PEFINEO grid FINITE ELEMENT CCNNECTION DATA. 

ELLIPS (APFGIL) 

GENERATES SUEREGICN CATA OVER REGULAR ANC DEFORMED ELLIPTICAL SHAPES. 
EXTRCT (REFINE) 

EXTRACTS REFIMeo GR(C FINITE ELEMENT CATA FROM REFINED NODAL DATA. 
EXTIO (BND) 

SAME AS EXTRCT, eUT FGR RCUNCARY FINITE CLEMENTS. 

GRDPHI (APFCIL) 

COMPUTES GRADIENT 8CUNCARY CCNCITICNS ALCNG AN SlRFOIl. SURFACE. 

MAXARA (ElKEY) 

DETERMINES MAXIMUM AREA OF A TRIANGULAR FINITE ELEMENT PAIR. 

POTENT (LINKK9)) 

SETS UP TGMPCRAPY STORAGE LCCATICNS IN IZ ARRAY FOR POTENTIAL FLOW. 
GUACR (PFFIME) 

PEPFCPMS EI-CUACRATIC T P ANS FORMAT ICN CF GU ADR ILATERAL 
SUBREGION CATA AND GENERATES REFINED GRID DATA. 

REFINE (APFCIL) 

PEPFGRMS GRID REFINEMENT CVES TWC-C I M ENS I CN AL SOLUTION DOMAIN. 

SFLBNC (REFINE) 

RENUMBERS THE SUEREGICN GENERATED BCUNCARY GRIDPOINTS. 

S EL INK (REFINE) 

FORMS SUEREGICN CONNECT ICM TABLE. 

SETOSN (REFINE) 

PERFORMS SUBREGION ELEMENT TO SUEREGICN NODE DATA TRANSFER. 

SHAPE (APFCIL) 

GENERATES SUBREGION COCRDINATES FOR SELECTED AIRFOIL SHAPE. 

SYMELM (POTENT) 

SETS UP CALL TC AIRFCIL CISCPETtZER ’ARFCIl.' ANC STORES OUTPUT 
■^0 BE USED FOR POTENTIAL SGLUTICN. 

TRIANG (REFINE) 

PERFORMS QUADRATIC TR ANS FOR M AT I CN CF TRIANGULAR SUBREGION DATA 
AMC GENERATES REFINED GRID DATA. 

XCOQRC (REFINE) 

TRANSFORMS SUBREGICN CCOPDINATES TO RECTANGULAR CARTESIAN. 


52 



TCPBULENT INTEC-R/L ECUNC/>RY LAYER SCLUT ION 


THE FCLLCWING SUEFCUTINES ^RE USEC TC COMPUTE UPPER AND LOWER 
AIRFOIL BCUNCARY LAYER INTEGRAL PARAMETERS. 


ACOE IFMAT) 

ARC (lAFLTI 
BLTINT ILINKUlOn 

COMPUTES TURBULENT INTEGRAL PARAMETERS TO BE USED FOR POTENTIAL SOLUTION. 
BOUND (BLTINT) 

CFCLY (ACCE.FAT) 

DRAG (ELTINT) 

FAT (ACOEiTURE) 

ENTRY - INFAT 
FMAT (INT8L) 

INIT (3CUND) 

INPUT (BLTINT) 

INSERT (TSLUl) 

INSTAE (TPNSIT) 

INTRL (TURBI 

entry - SETUPS 
LAMNAR (BCUNO) 

PERSON (INTEL) 

OUTPUT (TURR) 

PRNTER ( RUT IMT) , 

SCHORC (INPUT) 

SPUN (FMAT) 

SMOOTH ( INPUT ) 

ENTFY - SCERV 
STAG (BLTINT) 

TBLUl (BCUNC, FAT, INIT, I N PUT , L AMN A R , SMOOTH , TRC ALC ) 

TRCALC (TRNSIT) 

TRNSIT (BOUND) 

TUR5 (fiCUND) 

XSTEP (FAT, TURR) 

ENTRY - ISTEP 
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Program Variables 


I/SRR/SY 


THE FCLLCWING TS /> F/RTIflL LIST OF ENTRIES IN THE URRAY VECTOR 
WHICH CONTAINS INTEGER PARAMETERS THAT ARE USED TO 
CCNTPGL PROGRAM FLOW. NOT ALL ENTRIES ARE LISTED SINCE SOME 
GF THEM DC NCT PERTAIN TO THE FPGCLEM CLASS IN THIS DOCUMENT. 


lARPAY NAME OEFIMTICN NON-© DEF . 

ENTRY 

AR2 lAPPST - STARTING LGCA7ICN IN lARRAY CCATAINl.NG LENGTHS OF VECTORS 
TO EE ACCEC. 

S6 lOASE - BASE NC. FCR 12 ENTRIES. 200 

261 IBC - NUMBER OF BOUNCARY CCNDITICN TYPES. 

211 IPICPV - COCE TO ALLCCATE STORAGE FCR BI-DIRECTIONAL DERIVATIVE. 

GO IBL - 1 = POUMCARY LAVER PPCCPAM, 1 

0 ^ 2-C NAVIEP STOKES PROGRAM. 

203 IBOT - starting ELEMENT NC. FCR CEPUG PRINT IN DERVBL, GEOMFL AND STRF. 
255 ICIAGL - 1 = REMOVE OIAGCNALS FROM DISCRETIZATION PLOT. 

127 IDIFRT - NO. GF TIMES TC PRINT INTER. OUTPUT IN WLFLXS AND PRSGRD. 

491 ISXT - NO. OF EXTRA VECTORS TO BE ADDED TO IZ ARRAY. 

115 IFGPCE - 1 = REAC XCQR ANC YCQR IN 2E15.5 FORMAT IN SETUP. 

120 IFP - 1 = FROZEN CALCLLATICN IN THERMO. 

0 = NCN-FRCZEN CALC. IN THERMO. 

175 IFSL - WHEN ITKE = 1, the FCLLCWING CCNDITIONS APPLY, 

0 = INTEGRATE TKS AND DISSIPATION. 

1 = INTEGRATE TKc ANC USE MIXING LENGTH FOR LENGTH CONSTANT. 

2 = INTEGRATE TKE AND USE FREE SHEAR LAYER FOR DISS. LENGTH. 

03 IFSLT - TEMPORARY STCFAGE OF IFSL FOP. CALL TO GETFSL. 

123 IGAS - 1 = CALL CRHCGS. 

0 = CALL DRHCBL. 

Ill IMAX - NOCE AT which EOLINCARY LAYER THICKNESS IS FOUND. 

145 IM IN - SWITCH IN CKNLIN WHEN HHIN IS SET. 

158 INITCN - INITIALIZER IN CCMTES. 

28 IPASS - NO. OF CALLS TC DEPIV. 

250 IPFI - 1 = EXECUTE FCTENTIAL SOLLTICN. 

213 IFHIR 

271 IPRINT - 1 = PRINT CRCSS FLOW THICKNESSES IN BLTINT. 

105 IPTSPL - 1 = USE PATANXFR AND SPALDING'S FOR TAU WALL. 

0 ^ USE LLCWIEG - TILLMAN FCR TAU WALL. 

132 IPWPIT - COCE IN STRF FCR INTFRMEDIATE OUTPUT. 

IGO IRENC - END '^osmCN IN 'IZ* ARRAY. 

3 IROW - I OR 2 = INDEX FCR PRESENT OR PAST VALUE OF DEP. VAR. 1 

29 IRUN - PROBLEM NC. BEING RUN. ( USUALLY ONLY 1 ) 

259 ISIOE - NU.MBER OF SICES / SUPER ELEMENT. 

186 ITCA - UNIT NO. CN WHICH TO STORE IMT. PAR. DATA FOR PLOTTING. 

187 ITCB - UNIT NC . CN WHICH TC STORE 'PLOTS' DATA FOR PLOTTING. 

1?3 irOP - UPPER PRINT LIMIT CCLN'TER FOR ELEMENTS. 

97 ITKE - 0 = CC NCT INTEGRATE TKE - CISS. EONS. 

1 = USE TKE - DISS. TO CCMPLTE TU.RBULENT VISCOSITY. 

196 ITWALL - 1 = USE GUCY FOP TAU WALL. 

281 IVVY - DEPENDENT VAR. FCR WHICH TC CALL DPSISO. 

122 IWRIT - CEPUG PRINT FLAG IN DFCFBL, CCNTES AND WLFLXS. 
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494 IZPPC 

495 IZEROS 
92 IZSIZE 

492 IZSTRT 

253 JCOORC 

254 JSHP 


STARTING LCC4TICN GF iCCRESSES IN IZ ARRAY TO BE ZEROED OUT. 
NUMBER OF IZ ACDRESSES TC BE ZEROED OUT. 

VAXIMl.lf' Clf'ENSICN OF IZ VECTOR. 

STARTING LGCATICN IN IZ ARRAY VsHERE NEW VECTORS 
ARE TO BE ACOED. 

KEYS FGP INPUT COCRGINATE SYSTEM. 

1 = KARMAN-TPEFFTZ CLASS AIRFCIL. 

2 = CALL AFSFF ( USER SUPPLIFC SUBROUTINE TO SPECITY 

AIRFCI L CCCRCINATES. 


169 KCDC - I = RESET NLINE TO 50 ANC DUMP KODE TO 2. 

61 KOUMP - PRINT INPLT CAROS ANO CATA GENERATED IN 30INPT. 


4 KEVMTO - 


195 KFXBNC - 

25 KINO 

167 KNTPAS - 

6 KODG 

7 K0D5 

26 KCUNT - 
102 KGLT 

113 KFLVAR - 
C6 KENT 


5 2 K°OW 
136 KSAV 
272 KSKIP - 


INTEGPATICN TECFMCUE. 1 

1 = maximum absollie stability. 

2 = MAXIMUM FELATIVE STABILITY. 

3 = EULEP INTEGRATICN. 

FLUX PGUNCARY code USED IN WLFLXS. 

KIND CF ELEMENT. ( USED IN L7h). 4 

IF NO PRINT IN KNTPAS TIMES THRU OKNUIN, THEN FORCE PRINT. 

PRINT GECMETFY CUTPUT IF .NE. 0. 

PRINT INTER. CERIV OUTPUT KC05 TIMES. 

RUNNING CCUNT CF OUTPUT, t LIMITED BY LPRINT.) 

LOGICAL LNIT NUMBER CN WHICH TC STORE PLOT DATA. 

NO. OF VARIAELES to BE PLOTTED OR PUNCHED. 10 
PRINT CRTICN ( SET CURING EXECUTION.) 

0 = MO CALL TC FEOUTP ERCH QKNIMT. 

1 = CALL FEOUTP FROM OKMINT. 

MO. OF RCW5 IN D I SCR ET I Z AT I C N . 100 

PLOT TAPE NO. SAVED I.M OKMINT. 

N = PRINT EVERY N-TH INTEGRATION STEP IN BLTIMT. 


99 


50 LCCL 


47 LG 
213 LMLT 
179 LGC 
182 LOGS 
172 LCWO 
34 LPRINT 


NC. CF CCLUMNS IN CISCRETIZAT ION. 20 

IF .NE. 0 CN INPUT, THEN CNTPTS AND CNTNDS ARE 
TO BE READ IN . 

NO. CF CCLS. IN SCLUTICN FIELC. 

NQ. OF CCNTOLRS FOB WHICH TO COMPUTE MIX. LENGTH TURB. VISC. 
INTERVAL NC . FCUNC IN LOCK SUBROUTINE. 

PRINTOLT VAR. KCOE USEP IN STGUTl. 

USE LAMINAR VISCOSITY BELOW LOWO AND MLT FROM LOWD ON. 2 

LIMIT CN CLTFUT CCUNT. 100 


212 MLTRHS - NUMBER OF RIGHT HANC SIDES TO SOLVE FOR IN STRF. 1 

212 MLTRHS - COOE TO ALLCCATE CORE FOR MULTIPLE RIGHT HAND SIDES IN STRF. 

114 MSSC - CCNVERGEMCE SWITCH USED IN QKNUIN. 

103 MTM - NO. CF PASSES TG USE IN FECL'TP. (PROGRAM SET) 


256 NAFTP 
23 N8 
170 N8C 
121 NBORO 
69 N3SET 

33 NBUG 


NUMEER CF CATA POINTS SPECIFIEC ALONG ONE AIRFOIL SURFACE. 
NO. OF CHAP. IN EACH WCRO CF CUTPUT VAR. TITLE. 4 
MAX. NO. OF FCUNCAP.Y CONO . FOR ANY ONE OEP . VAR. 

NO. CF NCOES ARCUMO ECRCER CF C I SCRET I Z AT I ON . 

1 = STORE CEP. VARIABLE INTO YY ARRAY. 

0 = STORE YY ENTRY INTO DEP. VARIABLE VECTOR. 

1 THRU 5, PRINT DEBUG CLTPLT FRCM L7H AND POTENTIAL FLOW 
CATA GENERATOR. 
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22 NC 

125 NCALLS 

173 ^CC^'CC 

174 NCCi''TO 
59 NCFT/ie 

1 ND 
51 NCBL 
124 NOEfilV 
279 NfHM 
48 NDOF 
194 NOP 
168 NCPPES 
162 NOP\,SX 
14 NELEM 
89 ^E^'D 
31 NEQ 
43 NEC/5CC 

58 NECKNN 
45 NEXP 
107 NEIE2 


46 NF 

257 NFLUX 

258 NFX 
130 NGETh 

54 MHHALF 
53 NH2 
68 M 
63 NINO 
94 N!ZS 
65 NJ 
27 5 NL 
88 NLINE 
191 NM 


60 NMBGUT 
206 NM)L 


190 NMOUT 

16 NMGOE 
278 NNPT 
55 NODE 
150 NCDNO 
19 NOE 
148 NCNC 
142 NOUTPP 

283 NOl 

284 N02 
30 NP 

18 NPART 


- NC. OF CFAPAGTERS IN GUTFur FORMAT. 8 

- NO. OF ROUTINES TG CALL AT END OF OKMUIN. 10 

- NC. CF C7RC RF/C IN FOR COMCC TITLE PAGE. 

- NO. OF CARDS READ IN FOR TITLE INFORMATION. 

- NO. OF ENTRIES IN SPECIFIC HEAT TABLE. 1 

- IN'ITI ALIZATICN FARAPETER IN CfCFNS. 

- CKNUIN CODE TO DETERMINE MINIMLM STEP SIZE. 

-• 2 = CALL DERVEL . I 

- CORE ALLOCATION FOR DELTA STAR SUBROUTINES. 

- NC. OF DEGREES OF FREEDOM. I USED IN L 7H ) 1 

- SPACE ALLCCATICN IN INFINT VECTOR. 7 

- 1 = PRINT OUTPUT FROM PRSGRO BLT DO NOT USE IT IN SOLUTION. 

- NO. OF CFEX'S IN FFRIME TABLE. 4 

- number of ELEMENTS IN SCILTICN. 

- STAR-^ING LCC. PAR. IN FECIMN. 

- NO. CF VARIAELES TG BE SCLVED. 5 

- NO. OF ECNS. TO ADO AFTER TKE - DISS. STARTUP. 

- E.G. -2 = CELay INT. tkE ANC CISS. UNTIL C4EDSW IS SATISFIED. 

- NO. OF DEP. VAR. TG BE INTEG. IN QKNUIN. 1 


- NO. CF BCUNCAP NODES IN JBCC'NO VECTOR. 1 FROM L 7H ) 

- 0 = 00 NOT LSF MIXING LENGTH THEORY FOR DIFF. COEF. 

1 = USE MLT FOR SOLUTION OF DIFF. CQEF. 

2 = DELAY LSING MLT UNTIL FlE2S'rJ IS SATISFIED. 

- NO. OF 'MS' BYTE WORDS IN TITLE FOR EACH DEP. VAR. 4 

- NUMBER CF FLUX I A3) BOUNCARY CCMDITIONS. 

- NUMBER OF FIXEC NCOES. 

- COUNTER IN CPFCGS TO INIT. VARIABLES IN GAS. 1 

- NO. OF DECREASES IN STEP SIZE IN QKNUIN. 

- NC. OF INCREASES ’IN STEP SIZE IN CKNUIN. 

- STARTING LCC. IN CEP. VAR. MATRIX FOR THIS VARIABLE. 

- FEOIMN DISPLACEMENT CCUN'TEP. 

- NO. CF IZ FN'TRY POINTS THAT CAN BE STORED. 400 

- STARTING LCC. IN YY MATRIX FOR CEP. VAR. 

- MO. OF X, Z CGQPDINATES DESCRIBING LOWER SURFACE IN BLTINT. 


- LINE CCUNT CLTPUT CCNTPCL. 60 

- TYPE OF elements IN SOLUTICN. 3 

- 2 = LINE (CNE-CIMENS lOMAL) . 

- 3 = TRIANGLE (T kO -C I MENS ICNA L > . 

- NO. OF VARIABLES TO BE PRINTED. 30 

- ALLOW EXTRA STCRAGE IN [ 7.(711 ANC IZ(72) 

LENGTH OF IZ(71) = MAXIMLM ( NODE, NEO*NM*MMDL ) 
LENGTH OF IZ172) = MAXIMUM ( NODE, NEO*NMOL) 

- 3 = PRINT GL1PLT IN GECMETPY FCRM. I 

2 = PRINT OUTPUT IN NODE NO. SEQUENCE. 

- NUMBER CF NCCES IN SCLUTICN. 

- VARIABLE DIMENSICNING PARAMETER IN FEOIMN. 100 

- NODE NO. AT WHICH DPHOGS IS CCMPUTIMG. 

- NC. OF EGUAT 1CN BEING SOIVED FOR DEP. VAR. 'NP'. 

- NCN-CCNV^PGENCE CCCE IN GAS. 

- NO. OF SCALARS TO PRINT IN GLTPUT. 60 

- VECTCP FOP STCRAGE OF C/DX IN CPSISQ. 

- VFCTCR FOR. storage OF D/DY IN CFSISG. 

- CEP. VARTAELE BEING SOLVED AT THIS TIME-. 

- NO. OF PARiniCNS ( LSED IN ’ L7H' ) 2 
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198 NPGRCT - ST/!RTUP CnUMTER LSED IN PRSGRO. 4 

199 NPGBDV - ST/SRTUP CCUNTEP USEC IK PRSGRO. 4 

20 NPRNT - MO. OF PRIM PCSITICKS GK A LIKE CF OUTPUT. 132 

35 KPSICC - MO. CF SCUKC^oy NODES FOR VAR. ENTERING STRF. 

11 MPTOOF - NO. CF PCIMS CEG. CF FPEECON. ( JDCUND PAR.) I 

24 NPTEL - NO. OF PCIMS ELEMENT. (L7HJ 3 

153 KPUKCH - SET = 7 IF ELEMENTS AMD KOCES ARE TO DE PUNCHED IN DIMEN. 

161 NPVSX - NO. OF orpssLPES IN P VS X TABLE. 2 

49 KCEL - KO. CF PCIMS / ELEMENT. IL7H) 3 

21 MPOW - DEP. VAR. ANC CERIVATIVE ALTERNATOR IN QKNIMT. 2 
9 NPSTRT - LOGICAL TAPE KC. TO READ RESTART CONO. IN 'BDINPT*. 

£7 NRTAPE - LOGICAL TAPE KC. TC STCP6 RESTART CONO. IN 'LINK2' 

67 NS - GENERAL DUMMY PARAMETER. 

252 KSELEM - NO. OF SUPER ELEMENTS IN AIRFOI. 

146 NSFCBE - RESET CONDITICN FLAG IK 'FINOBE'. 

27 KSKIF - KO. OP EOUKCARY LQC. / DEP. VAR. NODE 

64 NSM - STOP program IF OLTPLT EXP. IS .GT. NSM. 10 

251 NSMOOE - NO. OF SUPER NODES IN AIRFOIL. 

121 NSPEC - KC. CF SPECIES IN SOLUTICN. 9 

273 NSTAG - 1 = INPUT GECMETPY AND PRESSURE DO NOT BEGIN AT STAGNATION 

FCIKT ANC EMC AT TRAILING EDGE. 

IBO NSTRT - starting LCCATICN IOF CEP. VAR. NCS . FOR SPECIES. 

276 KUS - NUMBER OF POINTS TC USE FOP COMPUTING DELTA STAR IN BLTINT. 
IC9 NTAFER - LOGICAL UNIT KC. CF RESTART TAPE. 

UO NTCNTS - STARTUP PARAMETER IN CONTES. 

62 KTITL - NO. CF TITLE CARDS TO EE READ IN AND PRINTED AT THE 
BEGINNING OF EACH OUTPUT SET. 10 

197 NTFPNT - 99999 = CC NCT PRINT INTEGRAL PARAMETERS IN TRBTHK. 

176 NTKS - NO. CF IKTEGPAL FARAMETEPS TO EE COMPUTED. 5 

ONLY 5 INT, parameters ARE COMPUTED, BUT OTHER PARAMETERS 
ARE SET, SUCF ANC TpETA REYNCLD’S MO., SHAPE FACTOR, ETC. 

IF = S9999, DC NOT PRINT INTEGRAL PARAMETERS BET'WEEN OUTPUTS. 

274 NU - NO. OF X. Z CCCPDINATSS CESCRIBING UPPER SURFACE IN BLTINT. 

232 NV - ARRAY FCR lEMCCPARY STORAGE OF CEP. VAR. IN DPSISQ. 

260 NVAP - NUMBER OF VARIABLES TO BE DISTRIBUTED OVER REFINED GRID. 


84 

NVARD 

- 

CCUNTER USEC 

IK STRF. 




85 

NVARl 

- 

COLNTEP LSED 

IN STRF. 




74 

NVH 

- 

CEP. VAR. NO 

. FCR ENTHALPY. 

DEF . 

= 

4 

7C 

NVP 

- 

DEP. VAR. KC 

. FCR FSI. 

CEF . 

= 

5 

71 

NVU 

- 

CEP . V AR . NO 

. FDR U1 VELOCITY. 

DEF. 

= 

1 

72 

NVV 

- 

DEP. VAR. NO 

. FOR U2 V ELEC ITY . 

DEF. 

= 

2 

73 

KIVW 

- 

DEP. VAR. NO 

. FOR U3 VELCCITY. 

DEF. 

= 

3 

90 

NYY 

- 

NO. OF TIME 

FERIGOS TO STORE YY. 

MUST 


2 

91 

NZZ 

- 

NO. CF TIME 

FFRICCS TO STCRE ZZ. 

MUST 


2 

193 

NM2 

- 

MMX'«2. LSEC 

FCR STORIMG FULL MATRICES. 




192 N2M - ,\M * 2. USEC FCR STORING SYMMETRIC MATRICES. 
166 N3CPNS - 1 = CALL PRSGRO FOP DPOX COMPLTATION. 
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R/lRR/iY 


T^E FCLLOWING IS A PARTIAL LIST OF ENTRIES IN THE RARRAY VECTOR 
WHICI- CCNTAINS REAL FAFAMETERS TFAT ARE USED TO 
CONTROL PROGRAM FLOW. NGT ALL ENTRIES ARE LISTED SINCE SOME 
GF TFEN CC NGT PERTAIN TO TFE PROBLEM CLASS IN THIS DOCUMENT. 


RARRAY NAPE 
ENTRY 


DEF IMTICN 

NQN-D DEF. 

156 

AINF 


REFERENCE SPFEC CF SCUNO. 


5 

AJ 

“ 

JOULES CONSTANT. 

778.28 

3 

ALC 

- 

CHARACTERISTIC ELEMENT SIZE. DEF. = 

MIN. SIDE 

2E 1 

ALPHA ~ 

- 

ANGLE OF attack. 


37 

ARNEW 

- 

NEW AREA CCMFLTATICN IN PRSGRO. 


86 

AVD 

- 

DAMPING FACrCR IN CFCFBL. 

25.3 

252 

BETA 

- 

KARMAN-TPEFFTZ camber AN8LE. 


209 

BLTh 

- 

BOUNDARY LAYER THICKNESS, DELTA. 


176 

CBTCK J 


SPECIFIC FEAT BRITISH TC NKS 

4.184 

365 

CD 

- 

TKE - CISS. CGEF. 

0.09 

36A 

CK 

- 

TKE - CISS. CGEF. 

1.0 

2 11 

CF0V2 

- 

SKIN FRICTICM 


18A 

cktke 

- 

TKE - DtSS. CGEF. 

0.09 

83 

CCPPX 


CCMPRESSICN FACTCR FOR CUTPUT COL. VECTOR 

INDICATES PERCENT OF XI AXIS TO BE USED TO SHORTEN 
SPACING INTERVALS. 

86 

COMFY 


CCMPRESSICN FACTOR FOR CLTPLT ROW VECTOR. 
SAME AS CCMPX, BUT FOR X2 AXIS. 


124 

COM 

- 

KARMANN'S CCNSTAMT LS60 IN MLT IN DFCF8L. 

.435 

70 

CCNFFC 

- 

IF .GT, 0.0, SET ALL RHO = CONRHO, 


62 

CCNV 

- 

CUTPUT SCALE FACTCR = 1 .0 / R£FL. 


77 

CGMl 


ALC / 1 RE*CFC INF*XHUINF ) 


78 

CCN2 

- 

CCM f TCFINF 


153 

CPA 

- 

SPECIFIC HEAT OF AIR. 

0.24 

159 

CPF 

- 

SPECIFIC FEAT CF HYDROGEN. 

3.445 

160 

CP INF 


SPECIFIC heat ccmfuteo in CPINIT. 


30 

CFOIMF 

- 

REFERENCE SPECIFIC HEAT. 

0.24 

153 

CVCF 


SPEC. FEAT CCNV. USED IN THERMO. 

4186.0 

148 

CVH 

- 

ENTHALPY CCNV. USED IN THERMC. 

1 .0 

151 

CVP 

- 

PRESSURE CCNV. USEC IN TFERM.O. 

.4725E-3 

152 

CVRHO 

- 

DENSITY CCNV. CSEC IN THERMC. 

16.02 

150 

CVT 

- 

TEMPERATURE CCNV. L SEO IN THERMO. 

1.0 

149 

CVU 

- 

VELOCITY CCNV. USEC IN THEPf-'O. 

0.3048 

132 

C ITKE 


TKE - DISS. CCEF. 

1.45 

133 

C2TKE 

- 

TKE - CISS . CCEF. 

0. 18 

120 

C4E0 

- 

MIXING LENGTH MULTIPLIER. 

0.0007 

143 

C4ECSW 

- 

• TKE ' DISS. STARTUP POSITION IN DFCFOL. 

30000.0 

13 

DELP 

- 

PERCENT INTERVAL FCR PRl.VTCLT. DEF. = 

2.0 

205 

OELSTR 

- 

cisplacement TUCKNESS. 


2C7 

CELT A3 

- 

ENERGY DISSIFATICN THICKNESS. 
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103 DEPLT - PERCEMT GF TG tq bF LSEO FOR PLOTTING STATIONS. 101.0 
165 CRTCDK - CFGPEES PA^K1^E TO GEGREES KELVIN 5.0 / 9.0 


175 EOTCKJ 
<50 EKMMF 
20<t ENER 
108 ENPULT 
14 EPS 
89 EPS INF 
95 EPTEST 
68 EP4M0 
368 EECF 
145 EIE2SW 


ENTHALPY BRITISH TG PKS 2.3244 

TKE NCN-G FACTOR - UINF*«2 
ENERGY FOR VELGCITV. 

GI PENS ICNAL I i ING FACTOR FOR ENERGY. 

ACCURACY TEST RET^EEN PRED I CT CP-CCR RECTOR FORMULAS. 0. 

GISSL. NCN-n FAC^GR - UINF'-«3 / ALC 

ZERO TEST FCF OISSIPATICN USEC IN DERVBL. 

MULTIPLIER FOR >HOCT IN PRSGRD. 1.0- 

SCALE FACTOR IN CISS. LENGTH . .435 

STATION AT LhICH TC CHANGE NE1E2 FPCM 2 TO 1- 30000.0 


1 FACT 
80 FACTH 
59 FACTPL 
79 FACT? 
370 FRCUCP 

163 FTTOCM 
162 FTTCIN 

164 FTTCMT 
189 FI 

329 FIO 


NCN-OIM. FACTCP. ( BL = ALC, 

1 . ) / ( CPC INF>»TQF INF) 

RHCINF « U INF A ALC 
l.C / FACIMU 

SCALE FACTCP USEC CN FROUOE NO 
FEET TO CENTINETEPS. 

FEET TO INCHES. 

FEET TC N ETEPS . 

Y-COOP.OINATE GF FI CUPVE. 


20NS = ALC / UINF ) 


. IN DERVBL. 

30.48 
12.0 
0 .3048 


LAST VALUE OF FI CURVE. 


31 G 

60 GAMraF 
142 Gl 
381 CIO 

187 G22 

188 G23 
140 G32 
14 I G33 


GPAVITI4TICN CONSTANT. 

FACTOR USEC IN GAS LAWS. 
Z-CCOROINATE CF GL CURVE. 

LAST value of Gl. 

VARIAELE GECNETPY FACTOR. 
variable GECPETRY scale FACTOR. 
VARIABLE factor. 

VARIABLE FACTOR. 


32.174 

1.4 


15 H 

121 HOOT 
88 HHMTN 
<57 HINF 

16 HMAX 

17 HM IN 
132 HRCCN 

122 HRSGCT 
45 HS 

7 HSIMT 
12 HT 
186 H21 
139 H3l 


CURRENT TRIAL STEP SIZE. 

Q MASS FLOW COMPUTED IN H2MIX. 

USEC IN CKNUIN FOR TIME STEP DETERMINATION. 

FEE. ENTHALPY. COMPUTED IN CRHCGS. 1.0 

MAX. STEP SIZE ALLOWED. DEE. = 0.02 * TO 

MIMNUM TNTEGPATICN STEP SIZE. 

CONSTANT LSEC IN H2MIX. 0.029126 

MIXING EFFICIENCY COMPUTED IN H2MIX. 

CU<JRENT STEP SIZE. 

START INTECRATIGN STEP SIZE AT THIS VALUE. 1.0 E-7 
OUTPUT VAR. FOR TIME STEF = HS « FACT / REFL 
GRID GPOWTF SCALE FACTOR. l.O 

GRID GPCWTH SCALE FACTCR. l.O 


2 CNE 
18 COTEPS 
375 GSHISC 
362 0S12 
361 OS 6 


FPCGRAM CCNSTANT. 

ACCURACY TEST PARAMETER IN 'QKMJIN'. 
I .0 / H21’»T2 

1.0 / FACTCPIAL (NMtl ) 

1.0 / FACTCRIAL(NM) 
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174 PDFTOC - POUNDS/FT**3 TC GRAMS/CM**3 0.01602 

170 PDFTCK - FCUA0S/FT*#3 TC KG/N>**3 16.02 

36 PF.CCIM - DIMFNSICNAL FRFSSUR6 = PEOGE * PRSCCN. 

39 PEDGE - NCA-CIM. PRESSURE AT PRESENT STATION. 

9 PIMP - FREESTREAf< PRESSURE. CEF. = 1ST VALUE IN P VS X TABLE. 

180 FI^SKGS - POUNDS / KG. 1.0 / 2.2 

19 PNTEFS - ACGURACY TEST PARAPETEP IN • QKNU IN • . 

99 PPRCCN - PRSCGN / ALC 

100 PPRIME - PRSSURE GRADIENT CCMPUTEO IN PRESSURE ROUTINE. 

67 PR - PRANDTL NUPGEP. l.O 

185 PRDIS - DISSIPATICN FRANOTL NUMBER. 1.3 

171 PPSCCN - RHCING ♦ UINF**2 / G 

181 PRTKE - TKE PRANDTL NUPRER. 1.0 

166 FSFTQA - POUNDS/FT+^2 TO PSIA 0.4725E-3 

169 PSFTOI - PCUN0E/FT*42 TC PCLNCS/1N**2 0.006924 

168 PSFTgn - P0UNnS/FT*^2 TG NE Vi TCNS/ P**2 47.38 

167 PSFTGT - PCLiN0S/FT*»2 TC TCFP 0.3591 

20 PTIM - PRINT TIPE PAPAMETER IN 'OKNUIN*. 

1C6 FTFL - FLCT point CEFIMTICN WFEn station PASSES PTPL. 


123 QR - CYNAMIC FPESSLP6 RATIO. . 1.0 

134 Q3MAX - MAXIPUP G CCNC. ECUND AT PRESENT STATION. 


179 RACCCN 
94 RAT02 
21 RE 
43 REEL 
47 REFLRE 
10 PHCINF 
157 RHOUIN 
8 RMINSl 
255 RNB 
199 PNULCC 
105 RGCST 
119 ROUALC 

116 PR 
2C8 RT 

32 RTCGNl 

33 RTCCN2 

34 RTCCNB 

55 RTCCN4 

56 RTCCN5 

57 RTCCN6 

117 RTOHMl 
28 RUN IV 


CCNVSRSICN FACTOR RACIANS TG DEGREES. 57.3 

RATIO OF Tva REMAINING GASES WHEN RUNNING H2 , 02 AND N2. 
REYNOLD'S NO. RHO JNF*U I N'F*ALC /XMUI NF 

REFERENCE LENGTH. UINE 1.0 

REYNOLD'S NO. RASED ON REEL.' 

FREESTPEAM CENSITY. 

RHOtNF ♦ UINE 

PRCGRAP GCNSTANT. -1.0 

P.EYNCLO'S NUPEEP FCP AIRFOIL. 

1.0 USED LOCAL VISCOSITY FCR VAN DRIEST DAMPING FACTOR. 
XMA / XPF - 1.0 

RHCINF * LINE * ALC**2 
CPF / CPA 

THETA REYNCLC'S NUPBER. 

2.') * G * AJ 

GAPPAF / 2 .0 

RUN IV / XMA 

RTCCN2 « XPACFG**2 

LINF**2 / (PKCM^CPCINF’STCFINF) 

2.0 •* RTCGN4 

RR * ( TCH/TCA - 1.0 ) 

UNIVERSAL GAS CONSTANT. 1545.33 


256 SANGLE - SWEEP ANGLE (CEGREES). 

129 SCT - CONSTANT SCHMIDT NUMBER. 

210 SH^PEF - SHAPE EACTCP. 

190 SLOPE - SLOPE OF VARIABLES CCMPLTED IN LOOK. 

6 SCUND - SPEED CF SCUNC FOR NODE BEING PROCESSED. 

Ill SPLIT - CUTOFF LSED IN H2MIX. 

44 SC2 - SORT! 2 .0 ) 

73 STLOCR - REF. CCN. TEPF. IN SUTHERLAND. 

74 STLDEX - EXPONENT USED IN SUTHERLAND. 


1.0 


0.02835 

204.0 

1.5 
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72 STLGTP - PEF. TENF. USED IN SUTHERLAND. 
71 STLDVR - VISCOSITY LSEO IN SUTHERLANC. 
50 SSINIT - FSIMIT / FACT 


A92.0 
. U63E-4 


374 TAREA - TOTAL CGMPLT47ICNAL AREA. 

372 TEAR - MASS WEIGFTEC AVFRAGE TEMPERATURE. 

35 TO - TOTAL SCLLTICN TINE ( CISTAACE J FROM TO. I .0 

22 TF - FINAL TIME {CISTftNCE), TF = TC + TO 

253 TEANG - TRAILING ECGf ANGLE. 

20C TFETA - MOMENTUM THICKNESS. 

196 TFETAP - SCALE EACTCP USED IN PRSGP.C. 

4 THK - OEF. NCN-CIM. THICKNESS CF ELEMENTS. 1.0 

254 TFKAF - AIRFOIL MAXIMUM THICKNESS. 

23 TIME - CURRENT TIME (CISTANCE). 

48 TIMESV - SAVED TIME LCCATICN FOR IMPLICIT INTEGRATION. 

24 TO - STARTING TIME (CISTANCE). 

146 TOA - AIR OEFEPENCE TEMP. FCR CCMPUTATICNS IN DIMEN. 533.0 

58 TCFINF - REFERENCE TEMPERATURE. 533.0 

147 TCH - H2 REF. TEMPERATURE FCR CGMPUTATIQNS IN OIHEN. 520.0 

40 TRATin - 1.0 ^ (GAMMAF-1.0) « XMACHS**2 / 2.0 

257 TRIFOP - 1.0 = CG NCT SET TRIP LOCAT ICM UNTIL RTHETA .GT. 200.0 

258 TRIPUP - X/C VALUE CF TRIP LOCATICN IN ELTINT. 

155 TSINF - static tcmpfratURE CCMP. IN CPINIT. 

26 TWOPt - PI « 2.0 


371 DEAR 
203 UGD 
63 UECCE 
27 UINF 
202 UWALL 


MASS WEIGFTEC AVERAGE VELOCITY. 
EOCE VELOCITY. 

UINF / UINFX LSEO IN BROSHW. 
FREESTREAM VELCCITY. 
velocity JLST CF WALL . 


104 VELCST 

177 VLBTCN 

178 VLBTQP 
IU2 V S TART 


UINF**2 / ( i.C 4 G + AJ # CPA * TCA ) 

VISCOSITY ERiriSF TO MKS . 1.488 

VISCOSITY ERITISH TO CGS . 14.88 

PERCENT OF TC AT WHICH TO START U2 COMP. IN CONTES. 


333 WSMAX - WIDTH CF OCMAIN AT INITIAL STATION. 


366 XEM 
125 XLAM 

37 XLE 
IG9 XMA 

61 XMACFO 
154 XMACHS 
373 XMDOTC 
66 XMF 
172 XMFACT 
110 XMH 
118 XMSOF 

38 XMUINF 
46 XMXOR 

193 XNWGEO 
9 8 X P R I M F 
52 XSCALE 
41 XSFFT 


EXRCNENT FACTCR USED IN VAN-CRIEST DAMP. FACTOR. 1.0 


CONSTANT USEC IN OFCFBL. 0.09 

LEWIS NUMRER. 1.0 

MOLELAR WEIGHT CF AIR. 28.97 

MACE NUMBER. 

LOCAL MACH NUMBER. 

AVERAGE MAS FLCW. 

MCLECULAR WEIGHT OF FLUID. 29.4 

UINF * SORT ( XMA / ( TOFIJVF * GAKMAF * G * RUN IV ) > 

MOLECULAR WEIGHT OF HYORCGEN. 2.016 

MASS DEFECT CCMFUTEO IN GMCCNC. 

FREESTREAM VI'COSITY. 


RATE OF CHANCE OF CPOX. (USED IN PRSGRD). 
l.C = RUNNING NVGECM. 

NGN-DIM. FRESUURE GRADIENT AT PRESENT STATION. 
XICCR SCALE FACTOR. 1.0 

SHIFT X-COCP.D INATE. 
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201 XTC - PRPSEM STiSTICN FOR INTEGRAL PARAMETER PRINT 

112 X3LAST - LAST TIME SIEF LSEC IN DERVCX. 


307 YLTKE 
330 YMULT 
198 YPLUS 
E3 YSCALE 
02 YSFFT 
76 YTT 


SCALE FACTCR IN DISS. LENGTH. , 4.5 

SCALE FACTCR FCP GRID MUITIFLIER. 

Y+ VALUE AT WHICH TO SWITCH FROM MLT TO TKE, 30.0 

X2C0P SCALE FACTCR. 1 .0 

SHIFT Y-CnCRD INATE . 

SCALE FACTCR FCR PLOTTING. 


302 ZVULT - SCALE FACTCR FCP GRID MULTIPLIER. 

29 ZT - DIRENSICNAL CLPPENT STEPIZE ICISTANCE). 
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IZ APP/Y 


THE FCLLCWING IS A LIST CF ENTRIES &T THE OEGIN’'IING OF THE 
IZ ARRAY WHICH CCMAIN ENTRY PCIMTS IN THE REMAINDER OF THE 
IZ ARRAY WHICH ARE THE STARTING LOCATION FOR THE VARIABLE 
length vectors LSEC in THE PRCGPAP. 


IZ 

ENTRY NAPE CEFINITICN 


1 

ICCL 

- 

2 

I ROW 

- 

3 

IFMTHO 

- 

A 

IT ITLE 

- 

5 

IIPINT 

- 

6 

IKBNQ 

- 

10 

1 I NCGL 

- 

11 

I INFOW 

- 

18 

ICPTA6 

- 

19 

1 TTAB 

“ 

20 

I ILSEC 

- 

21 

INPINT 

- 

25 

1 1 BMC 


26 

I INODE 

- 

27 

IJBCNin 

- 

28 

IKEYCL 

- 

29 

IK5Y0G 

- 

30 

IKFYRW 

- 

31 

IKIN 

- 

32 

INWN 

- 

33 

1 I NOEX 

- 

34 

1 INCRW 

- 

35 

INOCCL 

- 

36 

IlELS 

- 

37 

I lELEM 

- 

38 

I IBCRO 

- 

39 

IA2 llA 

- 

40 

1311 

- 

4 1 

182 11 

- 

42 

IP2 1 IS 

- 

43 

IC2CC 

- 

45 

ISMSTR 

- 

46 

IBCNST 

- 

4 7 

niF 

- 

48 

lYY 

- 

49 

IZZ 

- 

50 

IX 1P2 

- 

52 

lYCIM 

- 

S3 

IPCOL 

- 

54 

IPFOW 

- 

55 

lASUM 

- 

56 

lAHXLT 


57 

IWW 

- 


CISCRET IZAT ICN COLUMN LOCATICN'S. 

Cl SCRET I ZAT I CN RCW LCCATICNS. 

HEACINGS FCR CLTPUT VARIABLES. 

TITLE FCR start OF EACH CUTPUT PHASE. 

LIST OF OEP. VARIABLE NUPBEPS. 

NO. CF OGUNQARY NODES / OEP. VARIABLE. 

NC. GF NCDES FFP CCLUMN. 

NO. OF NODES PER ROW. 

SPECIFIC FEAT TABLE ENTRIES. 

TEMPERATURE TABLE ENTRIES. 

COUNTER USED IN FEOUTPUT . 

HEP. VAR. PC5ITICNS IN IIPINT VECTOR. 

RE-ORCER6C NCCES / DEP . VAR. TO ACCOUNT FOR BOUNDARY COND. 
ARRAY CF ELENENT CCNNECTICNS I NM/ CL EMENT ) . 

MODE SOLITICN ORDER USED IN BANCHC. 

CCLUMN KEYS FCP EANCHO. 

CIAGCNAL KEYS ECR E AHCKO . 

FOW KEYS FCR BNACHC.^ 

LIST OF CEF. ViR. NOS. FCR SPCCIES. 

TEMPORARY SrCFAGE FOP PE-CRCEPEC NODES. 

CRDER CF NCCES EY COLUMNS FROM LEFT TO RIGHT. 

CROER OF NCOES BY ROWS FPCM TCP TO BOTTOM. 

OUTPUT COLUMN PC'S IT ICN OF NODES BY ROWS. 

MO. OF ELE^E^TS CCNNECTEC TG NGCES. 

LIST OF ELEMENTS CCNNECTEO TO NODES. 

LIST OF ECRCER NCCES IN COUNTER-CLOCKWISE ORDER. 

A211A ANTI-STMMETRIC MATRIX. 

Ell MATRIX STCFAGF. 

R2LL MATRIX CCMFLTEO IN GECMFL. 

eziis matrix ccmputed in geomfl. 

B2C0 matrix storage. 

ENTRY POINTS IN IZ FOP, STANDARD MATRICES.. 

STORAGE FCR ECUNOARY CCNCITICNS. 

DIFFUSION COEFFICIENTS / OFF. VARIABLE. 

2 SPT.S OF VALUES / OEP, VARIABLE. 

2 SETS OF values CF CEPIV. / CEP. VARIABLE. 

ELEMENT lengths COMPUTED IN GECMFL. 

NON-CIMENSTCNAL transverse COORDINATES. 

COLUMN CnonOlNATES FOR CUTFUT FAGS. 

BOW CCCRCINATES FOR OUTPUT PAGE. 

SUM OF AREAS CF ELEMENTS AROUNC EACH MODE. 

MIXING LENGTHS . 

TE'-FCRARY STCFACE FCR CROSS VELOCITY. 
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58 IVV 

59 IVV 

fcl ITBfUK 

65 IU3P0S 

66 IU3VAL 

67 TVWALL 

68 IVViSTA 


T5MF0PARY STORAGE FOR TRANSVERSE VELOCITY. 

TEMPORARY STCRAGE FOR OOWMSTREAM VELOCITY. 

STORAGE FCF INTEGRAL PARAMETERS. 

DOWNSTREAM PCS. FOR TRANSVERSE COORDINATE CHANGE. 

SCALE FACTCR FOP TRANSVERSE COORDINATE CHANGE. 

WALL VALUE OF INJECTED TRANSVERSE VELOCITY. 

CCWNSTREAH STATION AT WHICH TO INSERT TRANSVERSE VELOCITY. 


71 lOL'Tl - 
TO 

76 IOUT6 - 

77 lAREA - 

78 ICP 

79 IH 
no IPSI 
82 10 
03 IQP 

84 IRHO 

85 ITEMP - 

86 IRHSF - 

88 ITK 

89 IXICCR - 

90 IX2C0R - 

91 IPRESS - 

92 I A ML 

95 INCRMY - 

96 INCRNZ - 

97 lUSQ 

98 lAVTHK - 

99 NJST 

I'lO IPRVAL - 

101 lOL 

102 lOPL 

103 IVEL 

104 IW 

105 IPRGRC - 
1G8 lYNGD - 

109 18112 - 

110 18113 - 

111 IPLOTS - 

112 IFLTYF - 
114 ISCHPT - 

117 INX 

118 INY 

119 I. MACH - 

120 lESTAT - 

121 ICALL - 

123 IGMULT - 

124 lOSAVE - 

125 INOUT - 

127 lICNCL - 

128 IICNDX - 

129 IWCOS - 

130 IVSIN - 

131 IIPAR - 


TEMPORARY STCRAGE. 

TERFCRARY STCRAGE. 

AREA OF ELEMENTS CCMPUTEO IM GEOMFL. 

MOCAL VALUES CE SPECIFIC FEAT. 

NOCAL VALUES CF ENTHALPY. 

TEMPCPARY STCPAGE'FOR VARIABLE TO BE SOLVED IMPLICITLY. 
temporary STCRAGE FOR OEFENOENT VARIABLE. 

TEM-ORARY STGRA8E. 

MOCAL VALUES CF DENSITY. 

NOCAL VALUES CF TEM-ERATLRE. 

RIGHT HANC SICE CF IMPLICIT EQUATION TO BE SOLVED. 

ELEMENT THICKNESS D I STR I BU T I CN . 

NOCAL VALUES CF TRANSVERSE COORDINATES. 

NOCAL VALUES CF NCPMAL CCORCINATES. 

NOCAL values CF -RESSURE. 

NCCAL VALUES CF LAMIN/ir VISCOSITY. 

NORMALIZED TRANSVERSE COCRDINATES. 

NCRNALIZEC VALUE CF CROSS COOR C UIATES . 

NOCAL VALLES CF PSI * 62115 * FSI. (UlUI) 

AVE. thickness of ELEMENTS AROUND EACH NODE. 

STORAGE FCP PHS IN IMPLICIT ECN . SOLVER. 

STORAGE FOR RESTARTING ’PRSGRO'. 

STCRAGE EOF CEPVCX ROUTING. 

STCRAGE FOR CERVOX FCUTINE. 

NOCAL VALUES C F U2 CCM-UTEO IN CONTES. 

NCOAL VALUES CF U3 . 

NOCAL VALUES CF DPOX. 

TRANSVERSE CCCPC. USED IN CONTES, OFCFBL, TRBTHK, WLFLXS, 
NATURAL CCCRCINATE DERIVATIVE CCMP. IN GEOMFL. 

NATURAL COCRCINATE CERIVATIVE COMP. IN GEOMFL. 

LIST OF VARIAELES TO BE -LCTTEC. 

TYPE OF PLCT TC BE GENERATED. 

NOCAL VALUES CF SCHMIDT NUMBERS. 

NUMBER OF CIVISICNS ALCNG XI CIRFCTION / SUPER ELEMENT. 
number of CIVISICNS ALCNG X2 CIRECTION / SUPER ELEMENT. 
NODAL VALUES CF HACH.NUMPEP. 

MOCAL VALUES CF STATIC ENTHAL-Y. , 

LIST CF LINK NCS. ANC ENTRIPS TO CALL AT END OF QKNUIN. 
LIST OF MULTI-LIERS FOP. CUT-LT VARIABLES. 

LIST OF VARIAELES TO BE PRINTED IN OUTPUT. 

TEM-ORARY’ STCPACE FOR OUT-UT VAR. AND SOURCE DATA. 

NO. CF MOOES / COL. USED IM CCNTCS, DFCFBL , TRBTHK, ETC. 
LIST CF NODES / CCL. USED IN CCMTES, DFCFBL, TRBTHK, ETC. 
COSINF OF ANGLE / COL. USED IM CONTES. 

SINE CF ANCLE / CCL. USED IN CCMTES. 

LIST OF -ARAMETEPS TO -PINT AT START OF OUTPUT. 


ETC. 
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132 IOP#R 

133 ISLTLC 

134 IPR 

135 r.VPfR 

136 !5PS 

137 ISKNFP 

138 I STM 

139 IX3ST 

140 I-VSO 

141 I0PX3 

142 rcpvx 

145 lAOIF 

146 IFXLT 

147 IXYCP 
197 lOOUTP 

201 NZCNF 

202 ISHAPE 

203 rCLZ 
2C6 M 
2C7 HYPE 

208 NCCL 

209 PEI EM 

210 NELK 

211 Q 

212 KC 

213 KF 

214 KTK 

215 IGCTYP 

216 OCUND 

217 ICUP 
2 18 BC 

219 MOCF 

220 KTFrX 

222 JELF 

223 KBAD 
22 5 A A 

226 BB 

227 SNOCE 

228 SGPIC 

229 CUN 

230 ELAK 

231 MAIL 


- title imfofmaticn for parameters at beginning cf output. 

- LIST OF CONSTiNTS USED IN SUTHLD. 

- NOCAL VALLES CF PRAMCTl NUNBER. 

- LIST CF ^ULT . FOR PARAM. AT START OF OUTPUT. 

- NOCAL VALUES CF TURB. VISC. CCN-UTED IN DFCFBL. 

- SKIN FRICTION CIST. BY CCLUMN. 

- STANTCN NC. CIST. PY CCLUMN. 

- LIST OF CCRNSTREAM STATIONS IN -RESSURE TABLE. 

- LIST CF CCVtN -PESSURES IN -RESSURE TABLE. 

- LIST CF Oa^MSTREAM StaTICNS IN PRESS. GRAO. TABLE. 

- LIST CF DOWNSTREAM -RESS. ERAO. IN DPDX TABLE. 

- NOPAL VALUES CF TURRULENT VISCCSITY. 

- MIXING LENGTH LSED IN DFCFBL. 

- AIPFCIL ECUNCiPY COORCINATESt CP'S, PHI AND NODE NUMBERS 

- MATRIX STORAGE FOR BNACHC SCLVER IN STRF. 

- super element numbers. 

- CROER OF GENERATEC ELEMENTS IN REFINE. 

- t(jpe of elements DESIRED / SUPER ELEMENT. 

- NUMBER OF ELEMENTS ALCNG CCORCIANTE 3. 

- TYPE OF SUPER ELEMENT 1 = TRIANGLE. 

2 = CUADP. ILATERAL . 

- NUMBER or CCLLMNS IN SUPER ELEMENT. 

- SUPER ELEMENT TC ELEMENT. 

- TYPE CF SUPER ELEMENT LINK. 

- SUPER NODE TC ELEMENT. 

- G KEY. 

- F KEY. 

- THK KEY. 

- RCUNCARY TYPE KEY . 

- BOLNCAPY CCEFF ICI ENT S. 

- CUM^'Y NCDE CCUNTER IN REFINE. 

- GENERAT<=C ECLNCARY CONDITICN ARRAY. 

- VECTOR OF ECUNCARY MOOES. 

- FIXEC CEF. VAR. NCCE KEYS. 

- VECCTR OF ROLNCARY ELEMENTS. 

- SUPER ELEMENT ECUNCARY KEYS. 

- NCN-PECTANGULf P C£RT-SIAN CCORC INATE5 DIRECTION 1. 

- NGN-RECTANGULAR CARTESIAN COORDINATES DIRECTION 2. 

- SUPER NCOE NUNEER VECTOR. 

- VARIABLE ARRAYS RETURNED FROM REFINE. 

- DUMMY STOP AGE liN REFINE. 

- SUPER NCCE CCNNECTICN TABLE. 

- GENERATED SUPER ELEMENT tq ELELENT DATA. 


201 - 240 ARE TEMPORARY LCCATICNS FCR PCTENTIAL AND DELTA * SOLUTIONS 
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Diagnostic Print 

Aside from the previously described standard problem oriented print, 
various forms of debug type print are accessible to aid in checking out 
program changes and size limitations. These print are accessed via keys 
specified in namelist NAMEOlt This section describes the more useful of 
these and illustrates the print to be expected for each. 


KDUMP As previously noted, the KDUMP flag 

presents a data reflection together with 
a print of the array filled by the data. 

In addition, it provides a print of the 
starting location of each variable in the 
II array as dimensioned in FEDIMN . This 
print is valuable for determining the 
array sizes when adding new vectors to 
the IZ array. A sample print of this 
output is illustrated in Figure 11. 

KODG This flag prints finite element inter- 

mediate data formed in GEOMFL . The 
print can become quite cumbersome if 
there are many elements in the solution, 
and the flags IBOT and ITOP are utilized 
to define the range of element numbers 
to be printed. 

K0D5 During the integration process the 

derivatives are formed by assembly over 
the finite elements. A print of the 
assembled and reduced vectors is obtained 
by setting K0D5 equal to the number of 
prints desired. The flag is decremented 
by one for each pass through the deriva- 
tive routine. Flags IBOT and ITOP also 
apply to this print. An example of this 
print is illustrated in Figure 13. 

IPWRIT Equation solving is performed in sub- 

routine STRF and this flag set to 1 causes 
the global system vectors and matrix to 
be printed during equation assembly and 
solution. Flags IBOT and ITOP indicate 
the span of elements over which print is 
desired. 

* KPNT must also be set equal to 1 in NAMEOl to obtain any debug print. 

KPNT will be automatically set equal to 0 in LINK2-6. 
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Figure 11 IZ ARRAY Vector Starting Locations in 
Hexadecimal) and Decimal Form 
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Figure 12 Finite Element Intermediate Print (6E0MFL) 
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PMULT _ KI-HAV MIH3ML , _ 

lo ?. t k.?*- ' i - 1 , ;ai6H-;u i.'joiooe oo o.o 

1) ' . CIF;-;UL( 2 ) . = ..-3 . . 0 !FMijL ( 3) = -- 2 . A 1 03*1 -0 1 . 

CNVlUU 1) = -A . 27"i.?--CA CN'jMUU ?.) = -I.OH^E-j 3 CNVMtjL ( 3) = ^.23 E -0 7 

avECi'f' AVEk'~' deUJSQ TKEOIS rjisdis" _ "tkepro 

0. l0Ct2r--0 3 C. ICCA2E-C3 0. 700233-02 0.150666-03 -0 . lOOOOt-O I -0.A3 1036-03 0. 22215E-07 

F211 »AT2IX ?']“ Cjr’n^SirM. .. \ Z . . - i 

L. . . I . CC :i:C6 OQ ... 2 . __-.1.COOOOE OC 3 .. -1 . JOJOOc 00 a l.OOOOOEOa 

5 t CCNTP rijijT tCN TO I.FS 

I 1.3 7 2756,0 0 2 . 1.3 6 5:026 . .. ...... . . 

Z e’..’ . .Z.. Z.ZZ-_Z ,Z1”ZZZZ.^ .Z.Z . ,Z - - - - 
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Figure 13 Debug Print From the derivative Routine 
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IWRIT 


The IWRIT key is used to flag debug 
print in three subroutines. In CONTES 
it is used to print column vectors of 
density-velocity data which is inte- 
grated to obtain U2- In DFCFBL , various 
turbulence parameters are printed de- 
pending upon the turbulence model used. 
Debug print from WLFLXS lists various 
intermediate data calculated in the 
boundary region of the discretization 
and used in wall shear stress calcula- 
tions. 
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appendix a 


Potential Flow 

Standard Check Case Data Deck 
Goradia NACA0015 (Mod.) 



6_ POTEriTIAL FLOW SOLUTION 


TZnPF 

8 . 

_.9FENAMP 

10 GNAMEOl 

11 NODE = 390. LCOL = 100. KROM = 100, NIZS = 250, 

12 flPVSX..? NPC =60. .. NEO = 2, 18L = 0, 

13 MVP = 5, NPRNT = 132, lOIFRT = 0, NOUTPR = 70, 

14 . MC = 7, NMOOUT = 1, KOOG = 0, IWRIT = 0, 

. .. 15_ -MRUG_=. 1,. .... MSPEC = 3, KOUMP = 1, KMTPAS =10, . . 

16 I ARRAY! 271 ) = l , 10 

17 1ARPAYI201) = 5, IARRAY(2051 = 100. IARRAYI200) = 5, 

__18„ I ARRAY! 2501 = 1, 57, 12, 0, 0, 0, 10, 54, 1, 4, 4, 1,.. . 

19 IARRAY(271! = 1, 10, 1. 

20 IARRAY12411 = 1, 10, 1, IARRAYI263I = 100, 

21 !ARRAY(276) = 2001, . IARRAY(279) = 200, 

22 GENO 

23 GNAHE02 

24 _„.T0FINF = 533.,__ UINF^r.,!*,. RHOINF = 1., PINF =1., 

2 5 ’ XMUINF = 533., 

26 XMUINF=533., i LC = .001, THK = .001, COMRHO = 1., 

27. THKAF,.= . .118. . ALPHA = 6., COM.PX ? 2., . COHPY =2., 

28 " ALPHA = 8., BETA = 0., THKAF = .15, RNB = 1.4E07, 

29 GEMD 

^ 31LINK4 TT.9 ~.’.V~~T . GE?^^^ G-GO f 6- PofENTl AL FLOW SOLUtIon' 

32NX _-i T. MO. Or OIV. PER SUPER ELEMENT NORMAL TO AIRFOIL 

_33 _ . . ,12*3 . T ... ... , __ 

34NY -1 T NO. OF OIV. PER SUPER ELEMENT TANGENT TO AIRFOIL 

35 64557675546 T 

_ _36LINK2 14 T... (OSCRT.Zl 

3700NE . . T ■ ENO CF LITERAL DATA 

38 . ... - . , 

39CaHTITLE T ..TITI.E. PRINTED 3EL0W COMOC, SYMBOL ,. . 

40 PERTURaATION POTENTIAL FLOW OVER' MACA 0015 IMOD.) AIRFOIL 

41 DONE 

42DESr.RIPT. ...20.4 _T ..OE SC R I PT I VE, T I TLE AT. BEGINNING .OF HEADER OUTPUT. 

43 .MACA 0015 (MOO.) AIRFOIL, 8 DEGREE ANCLE OF ATTACK. 

44 ... . 

4500NE _ T END CF LITERAL DATA . . 

460ESCRIPT 332 T lOPAR PARAMETER TITLES FOR HEAOFR OUTPUT, 

47 REFEREN'CE E'lGLISH-FT ENGLISH-IN M-K-S C-G-S 

_ 48 _L E NG T H ...... FEET . .. ....I N,...._...., ... ..M... «... ... «... .CM. «... 

VToone _ ”."7™ . ' .t" end of literal data 

50MPARA .-.I ...T MULTIPLIERS APPLIED TO HEADER OUTPUT. ILDC. IN RARRAY). 

51 542, 2=2 162 164 163 T 

52IONUMO . .-i .. .._ ...T LOCATION IN RARRAY OF SCALARS TO BE PRINTED IN HEADER 

53 _ _ .999, 5*200, 999, 200 4«43 T 

540ESCRIPT 203 f TITLES FOR OUTPUT OF DEPENDENT AMD PARAMETER VARIABLES. 

55PEP.TURBAT lOH PCTEMTIAL 

560CNE .t .. END . OF , LITERAL DATA 

STIOSAVe” -I T OEPEMDEMT VARIABLE AMD PARAMETER ARRAYS TO BE PRINTED 

58 5240 T 

59I0HULT -1 .T MULTIPLIERS APPLIED TO lOSAVE ARRAYS (LOC. IN ARRAY). 

”"''60 2 T 

61 

62IPINT -1 _T DFPfNOENT VARI ABLE .NUMBERS 

63 .. 5 T 

64KB.no .5 I T POTENTIAL FIXED AT THESE NODES 

65ADD _ .... 

66 I T 

67D0ME T ENO OF LITERAL DATA 

68 .. 

69VYY 

70 192*0. T 

71VYYEM0 ... 5 

72LINK3 4 T NON - 0 1 MEN.SIONALI ZE DATA (0K:EM) 

73LINK1 3 T FINITE ELEMENT MATRICIfiS tC.FOMFL) 

74C0M0C T PRINTS COMOC SYMBOL, 

75EXIT - T end CF JOB 
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APPENDIX B 
Boundary Layer 

Standard Check Case Data Deck 
Bradshaw Relaxing Flow 
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^BRADSHAW 

530BR 

6FEMAMC T READS NAMELIST DATA 

7 CNAMEOl 


a 


NODE = 60, 

NM 2* 

NEOKNN = 1, 

NMOUT 

= 3. 

9 

10 

11 

12 

CENO 

CNAME02 

NPRNT = 60, 
KNTPAS = 99, 

KDUMP = 1, 

lURIT = 10, 

KROW = 

28, 

13 


TOFINF = 533., 

UIWF = 1 10. . 

UINFX = 112.2, 

PINF = 

22 72. 

14 

15 


TO = 3.917 , 
REFL = 8.3333, 

TO = -5, 
KNULOC = 1.* 

OELP = 20. 

VST ART 

= 8. , 


16 CEND 

17FEDIMN T DIMENSIONS ARRAYS 

IBLINKl 1 T SETUP NODE - ELEMENT GENERATION * OUTPUT TITLE. 

19LINK2 16 ' T 

20VX2SCL 

21 0. 30 . 1 1.2 T 

22NDECRD 

23 1 28, 1 1. 0 T NODES 1-26 IN SOLUTION 

26ELEH 

2500NE T END OF LITERAL SEOUENCE 

26C0MTITLE T TITLE PRINTED BELOW COMOC SYMBOL 

276RAOSHAW SOLUTION - /DENT 2600, (REF. 3) 

28D0NE T END OF LITERAL DATA 

29C0M0C 1 PRINTS CCMQC SYMBOL. 

30DESCRIPT 206 T DESCRIPTIVE TITLE AT BEGINNING OF HEADER OUTPUT. 

31 

32BRADSHAW SOLUTION HLT 

3300NE T END OF LITERAL DATA 

36DESCRIPT 332 T I OPAR PARAMETER TITLES FOR HEADER OUTPUT. 


35 

REFERENCE 

E?^GL I sh-ft 

ENGLISH-IN 

K-K-S 


C-G-S 


36 

LENGTH 

. . . FEET.. . 

.»• .iw.*.. 

. . • . M 



• 

37 

VELOCITY 

..FT/SEC.. 




.. 


38 

DENSITY 

. LBM2ET3. . 


. .KG/M3. . . 


... G/ C C * • 


39 

TEMPERATURE 

.kANKlNt . ♦ 


..KELVIN.. 




60 

enthalpy 

. OTU/LBM. . 


..KJ/KO... 




61 

FROZ. SPEC. HEAT 

■BTU/LBM-K 


.KJ/KG-K. . 




62 

VISCOSITY 

.L6M/ET -S . 


.NT-S/M2.. 


..POISE.. 

* 

43 

LOCAL PRESSURE 

. . .psr., . . 

...PSl. ... 

. . N T / M2 . . . 


.. .TORR*. 

* 

44 

LOCAL SOLUTICN 

• NO. 

. .0!>0XI . . . 

..ENERGY.. 


• MIX. EFF 


65 

Xl/LREF 

.OXl/LHEF. 

. EPSl LON. . 

.OXIM/LREF 

REFL 

REYNOLDS 

NO 

66 DONE 

T END OF LITE 

•'AL DATA 





67MPARA -1 

T multipliers 

AfPLIc.D TO HEADER OUTPUT. (LCC 

. IN 

RARRAY) 


63 


5<2, 2*2 162 

17,6 163, 3*2 

164 163, 3*2 

170 

174, 


4S> 


• 3*2 165 2, 

2 -1 75 3*2, 

3*2 176 2, 

3*2 

177 178, 


50 


2 2 169 168 

167, 3*2 lOU 2, 

5*2 T 





SlIONUMO -1 

52 

53 
56 

55 

56 

57DESCRIPT 2U3 
58U1/UREF 
59D0NE 
60I05AVE -1 
61 

62I0MULT -1 
63 

661PINT -1 
65 

66K0NO 1 

67 DOT TOM 

68 LINK 3 6 

69LINK1 3 

70L1NK2 6 

71 VX3ST 

72 

73VPVSX 

76 

75L1NK2 20 

76LINK2 6 

77LIiVKCA(L -1 

76 

790KNI NT 
ilOEXl r 
81CASE END 


T LOCATION IN RARRAY OF SCALARS TO BE PRINTED IN HEADER 
999, 5*200, 959, 200 4*6.',, 200 27 200 2*27 , 200 
10 200 2*10, 200 58 200 58 200, 200 97 200 97 

2U0, 200 30 2C0 30 200, 200 33 200 2*38, 

999, 39 6*36, 200 156 98 135 122, 

11 12 16 05 67 T 

T IFMTHO TITLES FOR OUTPUT OLPENOtNT VARlAbl.ES. 

U2/UP.EF Ul PRIME EEL. MU/MUREF 

T END OF LITFRAL SEQUENCE 

7 OEPEMDENT variable AND PARAKLTEr. ARRAYS TO BE PRINTED 
1268 2268 1269 1267 T 

T MULTIPLIERS APPLIfO TO lOSAvE ARRAYS I LOG . IN ARRAY). 
6*2 T 

T UEPENOENT VARIABLE NUMBERS 
1 2 T 

T Ul FIXED AT THESE NODES 
DONE 

T NON - OIHENSIONALI ZE DATA I 01 MEN) 

T FINITE ELEMENT MATKICIES (OEOHED 
T CONTINUITY EOOATIOM SOLVER ICuNTES) 

0.08333333, T XL C CURD 1 NAT ES FOR PRESSURE TABLE. 

67.0 53.0 59.0 65.0 200. u T 

T PRESSURE DATA FOR PRESS. TABLE (CP) 

1.04 1-02 1.00 1.00 1.00 T 

T GENERATES Ul INITIAL PRuFILL FOR BRADSHAW CASE 
T PRINTS GLNEPATED NODE MAP 

T PLACES CALLS 70 Ll.NKI J AT END . OF OKfAJIN. 

24, 15, 23, 56, 215, T 

T INITIATES I NF EGRAT ION RETURNS CUNT RUL TU bOINRT AT TF 
1 END OF JOB 
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APPENDIX e 


Boundary Layer and Wake Flow 
Standard Check Case Data Deck 
Joukowski , 12% Thick, 6° Angle of Attack 



I 


. 2 
3 







42DBL 






5FENAME 

T RFAOS NAHEtrST DATA 



_ 6 .. f.NAMEOl. 

. 


. 


^ 

7 . 

NODE = 55, 

KROH = 55, 

LCOL = 2, 

NI ZS = 

200, 

8 

NEO = 4f 

N5QKNM 4, 

NE1E2 = 1, 

•lEQADO 

- -3, . 

9 

N0ERIV.= 2. 

NH _ . =2, 

ITKE = 0, 

NTKS 

= 10, 

10 

NPVSX = 70. 

NBC = 2, 

NCUTPR=80*, 

KOUMP 

= 1 , 

11 

NMOUT 3, 

KNTPAS=49, 

IFSL = 0, 

LG 

= 2, 

.. 12__ 

NSCY .= .1, 

IPTSPL = 1, . 

lARRAY (206) = 

20, NC . 

.= .10, - 

13 

KPNT = 1, 

I WRIT = 1, 

ITDA = 0, 

. ITDB = 

0, 

. 14 . 

as f.ENO 

NU2POS=20» 

NU3P0S=20, 

NPPNT=60, 



16 CNAME02 

17 

TOFINF = 533.. 

PINF = 2116.8, 

UINF = 40., 

XMUINF^ 

=.1238E-4, 

18 

REFL = 1., 

. AOUCT =_ 1.. 

RARRAY(391) = 

.88.RARRAYI 392)=1.6667, 

19 

PRTKE=1. , 

PROIS = 1.3, 

CITKE = 1.45 

, C2TKE ■ 

= - 18, 

20 

CKTKE = .09, 

CO = .09, 

YLTKE = .435 

, ESCF = 

1,, CK=1., 

21 

C4EDSW „ .999, 

E1E2SW = 5. , 

YPLUS = 8., 

VSTART 

= 10. t 

22 

TO = 499i 

TO = .01, 

OELP = 20., 

HMAX = 

10000., 

23 

24 

HSINIT = l.OE- 

07, 




25FE01HN 

. . T DIMENSIONS ARRAYS 





26 ,, 

27LINK2 . T GcMERATe GRID FOR VISCOUS SOLUTION - 2D lOSCRTZ) 

28VX2SCL T INTERPRET SCALE FACTORS AND LIMITS. 

29 .97, 15 l.O 0.840, I 1,0 .1.0, 20 1.06 1.13 T 

. 30NQECR0 .. .T ...SPEXIFY NO. OF NODES. TO SE OENCRATcD .IN EACH OIR. . 

31 . 1 37, I 1, 0 T 

32ELEM . T CONSTRUCT F. E. CONNECTION TAOLE. 

3300MF T END OF LITERAL SEOUENCF 

34CNTPTS -1 T NUMBER OF NODES IN EACH SUPER ELEMENT 

35 21 16 T 

. 36CNTN0S . .T. .RENUMBERS THE MOOES TO l.MCREASE FROM..BNO. OUT. 

) 37 21*T 1 17, 16'»I-1 16 T 


38C0MTI TLE 

T 

T ITLE PRINTED 

BELOW COHOC SYMBOL 




39 AIRFOIL trailing EDGE AND WAKE 

SOLUTION 





4000NE 

T 

END OF literal 

DATA 





410ESCRIPT 204 

r 

DESCRIPTIVE title AT BEGINNING OF HEADER 

OUTPUT. 

42 . 








43 JOUKOWSKI 

T/C = 

.lie, 6 OEG. ANGLE OF ATTACK 

, BND.-HAKE S 

OLN. 


4400NE 

T 

ENO CF LITERAL 

data 





450eSCRiPT 332 

, T 

JOPAP PAP.AMET 

FR TITLES FUR 

HEADER OUTPUT 

• 


46 REFERENCE 


ENGLISH- FT 

EMCLISH-I N 


M-K-S 


C-G-S 

47 LENGTH 

48 V^lCCITY 


• . • FEFT . . . 

IN.... 


• • . M . • . • . 


a»aaC^*.*a 

.a aCM/Sa,* 


. .FT/SEC .. 


..M/S 


49 DENSITY 


,Lr>M/FT3. . 



•KG/ M3. , . 


a. aG/COa a a* 

50 TfiMPGRATURE 


.RANKINE.. 



.KELVIN. . 



51 .ENTIIALwy 


. BTU/L3M. . 



a K J / KG > • a 



52 FROZ- SPEC. HEAT 


.BTU/LHM-R 



KJ/KG-K a . 



53 VISCOSITY 


.LB.M/FT-S. 



NT-S/M2.. 


..POISE... 

54 LOCAL PRESSURE 


. . .PSF. . . . 

• * * PS I « • • • 


.NT/H2. a . 


a , alnPRa a a 

55 LOCAL SOLUTION 


.MACH. NO. 

. -DPDXl,. . 


.energy . . 


.MIX. EFF, 

56 NWGEOH M'S 


-. ,H21.,.. 

. . .G22. ... 


a.G23a«.. 


.,..Fl 

57 XI/LREF 


.OXl/LREF. 

.EPS I LGM. . 


DXIM/LRFF 

REFL 

PFTiNOLOS MO 

58nONF 

T 

CND OF LITERAL 

DAT 
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115 . _ . _ . _ . 

116RESTART 9 1 T RESTART AT TRAILING EDGE (PARABOLIC MAV. STOKES) 

U7 

118NAMELIST .. . . 

119 CNAHEOl 

120 NEOKNiN = A, MEOADO = 0, NOBL = 0, ME1E2 "0, 

121 . _ KDUMP = _ , 

122 SEND 

123 GNAHE02 

_12A _ _ HMIN = 0., ^ _ 

125 SEND , . . 

126IARP.AY . IA5 .0 T RESET IMIN TO ZERO 

. 127RARRAY __13 0.1 -3 _T RESET PRINT INTERVAL (NON-0) _ _ _ 

128RARRAY 15 7 -IT RESET STEP SIZE (NON-DJ 

129RARRAY 22 l.OOl -I T RESET FINAL STATION (NON-D) 

„130ELEMENTS _t1 T ..RESET finite ELEMENT COMMECTION TABLE 

131. 1 2 2 3 3 A 4 5 5 6 6 7 7 8 8 9 9 10 10 11 11 12 12 

132 _ . 13 13 14 14 15 15 16 16 17 16 18 18 19 19 20 20 

.133 21 21 22 22 23 23 24 24 25 25 26 26 27.27. 

134 28 28 29 29 30 30 31 31 32 32 33 33 34 34 

135 35 35 36 36 37 T REN. 17 18, 16 18 

. .1360ESCRTPT. 203 .T. .TITLES FOR OUTPUT OF DEPENDENT AND PARAMETER, VAR I ABLFS,..„_- 

137U1 / UREF U2 / UREF TKE / TKEREF OISS / DISSREF MU / MUREF 

138U1 PRIME U2 PRIME TKE PRIME DISS PRIME DPDX 

_ 139D0NE T END OF LITERAL DATA.. 

140I0SAVE -1 T DEPENDENT VARIABLE AND PARAMETER ARRAYS TO BE PRINTED 

141 1248 2248 5248 6248 1247 1249 2249 5249 6249 305, T 

_142I.0MULI .. _-l T_ .MULTIPL lERS. APPLIED .TO lOSAVE ARRAYS ILOC. IN ARRAY). .... 

143 ■ 10*2 T 

144K8N0, .1 T U1 FIXED AT THESE NODES 

145D0NE _ T END OF LITERAL DATA 

146 1 BORO T REORDERS BOUNDARY ELEMENT NODE PAIRS 

147RIGHT DONE 

148K8Mn_. 2 _1 T OU/DX BOUNDARY COMOITIC.N ON U2 , 

149T0P 6 d 0 0. 2 -4.321 2 T ' 

15080TT0M , 0 ,.,„,0 ,. 0 . 0. 2 4.321 .. 2 . T 

15100NE T END OF LITERAL DATA 

152KBNO 5 T TKE FIXED AT THESE NODES 

15300NE_ T END OF LITERAL DATA . .. 

154K8NO 6 T DISSIPATION FIXED AT THESE NODES 

15500NE T END OF LITERAL DATA 

156LINKCALL_ -I __T_ ...P L A C E S CALLS TO LINKI J. AT END OF QKNUIN. 1 

157 5 1, 1 4. 2 15, 5 6, T 

15830PNS 

.1590KNINT... ... T INITIATES INTEGRATION, RETURNS CONTROL TO BDINPT .AT TF 

160 

161RESTART 0 1 T RESTART TO INCREASE STEP SIZE INCREMENT 

162 . . _ ... 

163NAMELIST 
164 SNAMEOl 

. 165 _. ... .NOEL .= P, . ... 

166 SEND 

167 CNAME02 

.. 168.. . , HMIN. = ,.0.0, . . . . 

169 SEND 

170IARRAY 145 0 T RESET IMIN IN NOMEOl 

171RARRAY . 115 .7.-1 T RESET HMIN NON-DI MENS tONAL I ZED BY FACT . 

172RARRAY 20 23 T RESET NEXT PRINT IIME 

173RARRAY 1 35 0, 1 T RESET PRINT INTERVAL INC-N-D) 

174RARRAY 22 1.25 -1 T RESET FINAL STATION -TF- (MON-U) . 

175RARRAY 35 0,25 -3 T RESFT SOLUTION INTERVAL -TD- INON-0) 

176SAVETAPE 9 4 

1770KNINT . T INITIATES I MTEGRAT ION, RETURNS CONTROL TO BDINPT AT TF 

178EXIT T END OF JOB 
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